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Three-dimensional simulations for the merger of binary neutron stars are performed in the frame- 
work of full general relativity. We pay particular attention to the black hole formation case and 
to the resulting mass of the surrounding disk for exploring possibility for formation of the central 
engine of short-duration gamma-ray bursts (SGRBs). Hybrid equations of state are adopted mim- 
icking realistic, stiff nuclear equations of state (EOSs), for which the maximum allowed gravitational 
mass of cold and spherical neutron stars, Msph, is larger than 2Mq. Such stiff EOSs are adopted 
motivated by the recent possible discovery of a heavy neutron star of mass ~ 2.1 ± 0.2Mq. For 
the simulations, we focus on binary neutron stars of the ADM mass M ^ 2.6Mq. For an ADM 
mass larger than the threshold mass Mthr, the merger results in prompt formation of a black hole 
irrespective of the mass ratio Qai with 0.65 ^ Qm < 1- The value of Mthr depends on the EOSs 
and is approximately written as 1.3-1.35A/sph for the chosen EOSs. For the black hole formation 
case, we evolve the spacetime using a black hole excision technique and determine the mass of a 
quasistationary disk surrounding the black hole. The disk mass steeply increases with decreasing 
the value of Qm for given ADM mass and EOS. This suggests that a merger with small value of 
Qm is a candidate for producing central engine of SGRBs. For M < Mthr, the outcome is a hyper- 
massive neutron star of a large ellipticity. Because of the nonaxisymmetry, angular momentum is 
transported outward. If the hypermassive neutron star collapses to a black hole after the longterm 
angular momentum transport, the disk mass may be ^ O.OIM0 irrespective of Qm- Gravitational 
waves are computed in terms of a gauge-invariant wave extraction technique. In the formation of the 
hypermassive neutron star, quasiperiodic gravitational waves of frequency between 3 and 3.5 kHz 
are emitted irrespective of EOSs. The effective amplitude of gravitational waves can be ^ 5 x 10~^^ 
at a distance of 50 Mpc, and hence, it may be detected by advanced laser-interferometers. For the 
black hole formation case, the black hole excision technique enables a longterm computation and 
extraction of ring-down gravitational waves associated with a black hole quasinormal mode. It is 
found that the frequency and amplitude are « 6.5-7 kHz and ~ lO"'^'^ at a distance of 50 Mpc for 
the binary of mass M ~ 2.7- 2.9Mq. 

PACS numbers: 04.25.Dm, 04.30.-w, 04.40.Dg 



I. INTRODUCTION 

Binary neutron stars in close orbits are strong emit- 
ters of gravitational waves. A scenario based mainly on 
a recent discovery of the binary system PSRJ0737-3039 
suggests that the detection rate of gravitational waves 
by the advanced laser-interferometric gravitational wave 
observatory (LIGO) will be ~ 40-600 yr'^ 0; i.e., one 
event per year is expected within the distance of « 35-90 
Mpc since the advanced LIGO is able to detect gravita- 
tional waves from coalescing binary neutron stars within 
the distance of about 300 Mpc. This indicates that a coa- 
lescing and merging binary neutron star is one of the most 
promising sources for kilometer-size laser-interferometric 
detectors 

Merger of binary neutron stars has been also proposed 
for many years as a possible formation scenario for a 
central engine of short-hard gamma-ray bursts (SGRBs). 
Associations between SGRBs and elliptical galaxies re- 
ported recently [7] make it unlikely that SGRBs are re- 
lated to supernova stellar core collapse since elliptical 
galaxies have not produced massive stars in the past 

10^° yrs. In addition, recent observations of the af- 
terglow of the SGRB 050709 rule out the presence of a 



supernova light curve and point to a binary compact ob- 
ject merger as the most likely central engine 8, 9j]. The 
merger of compact-object binaries (binary neutron stars 
or black hole-neutron star binaries) is thus the most fa- 
vored hypothesis for explaining SGRBs. According to a 
standard scenario based on the merger hypothesis, after 
the merger, a stellar-mass black hole is formed with a sur- 
rounding accretion torus of mass > O.OIM0 Energy 
extracted from this system by either magnetohydrody- 
namic processes or neutrino radiation powers the fireball 
and high-Lorentz factor jets for SGRBs, for which the 
typical burst energy is between lO'** and lO'*^ ergs (after 
a correction of a beaming factor) 

Hydrodynamic simulations in the framework of full 
general relativity provide the best approach for studying 
the merger of binary neutron stars. Over the last several 
years, numerical methods for solving coupled equations 
of the Einstein and hydrodynamic e quat ions have been 
developed (e.g., jTHIICs, 14, ll El 111 El El IS 113) 
and now such simulations are feasible with an accuracy 
high enough for yielding scientific results (e.g., [rRl2l| ). 

For many years, the fully general relativistic (GR) sim- 
ulations for the merger of binary neutron stars had been 
performed adopting an ideal equation of state (EOS) 
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P = {V—\)pe where P, p, e, and F are pressure, rest-mass 
density, specific internal energy, and adiabatic constant 
[T^ IT3I E, E|- However, simulations with realistic 
EOSs are necessary for quantitative understanding for 
the merger. From this motivation, in a previous paper, 
we performed simulations taking into account realistic 
nuclear EOSs f2ll| and clarified that properties of the out- 
comes and gravitational waveforms are significantly dif- 
ferent from those obtained in the ideal EOS with F = 2. 
In particular, we found that a hypermassive neutron star 
|22ll23l | of elliptical shape is formed for a reasonable mass 
range of the merger progenitor because of the stiffness 
of the realistic EOSs. As a consequence, hypermassive 
neutron stars can be strong emitters of high-frequency, 
quasiperiodic gravitational waves, which may be detected 
by the advanced laser-interferometric gravitational wave 
detectors [2^ . 

In this paper, we extend the previous work from the 
following motivation. One is based on the recent pos- 
sible discovery of a high-mass neutron star with mass 
2.1 ± 0.2Mq (one a error) '25^. This measurement indi- 
cates that the maximum mass of spherical neutron stars, 
Msph, may be larger than ~ 2Mp,, and hence, very stiff 
EOSs such as those proposed in j2a, |23 are favored. In 
particular, if Mgph is larger than ~ 2.IM0, the nuclear 
EOS is so stiff that most of the EOSs proposed so far 
will be rejected, resulting in that restricted EOSs such as 
those derived by Akmal, Pandharipande, and Ravenhall 
(hereafter APR) '26] survive. However, fully GR simula- 
tions for the merger with such stiff EOSs have not been 
performed to this time. This paper is mainly devoted to 
presenting new numerical results obtained with the APR 
EOS. 

The other motivation is to determine the final outcome 
for the case that a black hole is formed promptly after 
the merger. Clarifying the final state of the black hole 
system, in particular the surrounding disk mass, is an 
important subject in exploring whether a remnant of the 
merger can be a central engine of SGRBs. The previous 
simulations [l^ I2H were performed only for a short time 
after the formation of black hole because of the so-called 
grid stretching around the black hole horizon. Under 
this situation, the geometry is so steep near the horizon 
that the accuracy in the numerical computation is lost. 
The popular approach for overcoming this difficulty is 
to adopt excision techniques 28,'2^|33- In the present 
simulation, the growth of the black hole is followed until 
the system approximately reaches a relaxed state by em- 
ploying a simple excision technique [2^ [s^l • With this 
technique, the disk mass surrounding a black hole can be 
determined. Simulations are performed for a wide vari- 
ety of mass ratio for the black hole formation case as well 
as in two EOSs. We show that the disk mass surround- 
ing a black hole depends strongly on the mass ratio of 
the binary. In addition, the excision technique enables 
a simulation long enough to extract gravitational waves 
emitted after the formation of black holes. We show for 
the first time that gravitational waves are determined by 



a quasinormal mode of the formed black hole. 

To determine the disk mass around a formed black 
hole, Oechslin and Janka recently performed a series of 
interesting simulations for merger of binary neutron stars 
im using a realistic EOS (the so-called Shen's EOS) [H, 
which is also so stiff that the maximum mass of spheri- 
cal neutron stars is « 2.2Mq. They employ an approxi- 
mate formulation of general relativity (the conformal flat- 
ness approximation for the spatial three-geometry). In 
this approximation, gravitational radiation is neglected. 
They follow the merger process only in the formation of 
hypermassive neutron stars. Assuming that the hyper- 
massive neutron stars eventually collapse to a black hole 
by a dissipation mechanism, they estimate the disk mass 
(i.e., they compute the rest-mass of fluid elements with 
the sufficiently large specific angular momentum which 
are expected to escape falling into the black hole) . They 
conclude that the disk mass around the formed black 
hole will be always larger than 0.02A/o because the an- 
gular momentum transport during the merger process 
efficiently works. As shown in Sec. IV, the angular 
momentum transport indeed plays an important role in 
the formation of disks around the hypermassive neutron 
stars, and hence, our numerical results agree partly with 
theirs. However, for the merger of high-mass binaries, 
our results do not agree with theirs. 

In their case, the merger results in a hypermassive neu- 
tron star even for the total gravitational mass « ZMq. 
This result disagrees with ours. The reason may be partly 
due to different choice of the EOS from ours, but is also 
likely due to the fact that they do not take into account 
radiation reaction of grav itational waves. As shown in 
Sec. as well as in |2lj . angular momentum is dissi- 
pated by a factor of ~ 15% in the first 3 ms from the last 
one orbit. This significant dissipation induces prompt 
collapse to a black hole for the high-mass case in our re- 
sults. In this C0jSG, Oil black hole is formed in much shorter 
than 1 ms after the onset of the merger, and hence, the 
angular momentum transport does not work efficiently 
before the black hole formation, unless the mass differ- 
ence of two stars is significant. We find that the disk 
mass is much smaller than O.OIMq in the prompt for- 
mation of a black hole for the nearly equal-mass case, as 
shown in Sec. IV. 

The paper is organized as follows. In Sec. II A- 
D, basic equations, gauge conditions, excision scheme, 
methods for extracting gravitational waves, and quan- 
tities used in the analysis for numerical results are re- 
viewed. Then, the hybrid EOSs adopted in this paper 
are described in Sec. II E. In Sec. HI, initial condi- 
tions and setting for simulations are described. In Sec. 
IV, numerical results are shown, paying attention to the 
merger process, the formed outcome, and the disk mass 
surrounding a black hole. Implication of our results to 
formation of a central engine of SGRBs is also discussed 
in IIV Dl In Sec. V, gravitational waveforms are pre- 
sented. Section VI is devoted to a summary. Through- 
out this paper, we adopt the geometrical units in which 
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G = c = 1 where G and c are the gravitational constant 
and the speed of light. Latin and Greek indices denote 
spatial components {x, y, z) and space-time components 
{t^x,y,z), respectively: r 



denotes the Kronecker delta. 
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II. FORMULATION 
A. Summary of formulation 

Our formulation and numerical scheme for fully GR 
simulations in three spatial dimensions are the same as 
in 0, 1^ , to which the reader may refer for details of 
basic equations. 

The fundamental variables for geometry are a: lapse 
function, (3^: shift vector, 7^: metric in three- 
dimensional spatial hypersurface, and Kij-. extrinsic cur- 
vature. In addition, we define the conformal factor 
tp = e'^ = 7"'^/"'^^, conformal three-metric 7^ = e~'''^7y , 
three auxiliary functions Fi = S^'^djjik, the trace of 
the extrinsic curvature K, and a tracefree part of the 
extrinsic curvature Aij = e~^'^{Kij — ^ijK/?,). Here, 
7 = det(7„). 

The Einstein evolution equations are solved using a 
version of the BSSN formalism following previous papers 
We evolve 7^, 0, Fi, Aij, and K using 
an unconstrained free evolution code. The latest version 
of our formulation and numerical method is described in 
[l7| . The point worthy to note is that the equation for 
is written to a conservative form similar to the continuity 
equation, and solving this improves the accuracy of the 
conservation of the ADM (Arnowitt-Deser-Misner) mass 
and the total angular momentum significantly. 

The fundamental variables for the hydrodynamics are 
p: rest-mass density, e: specific internal energy, P: pres- 
sure, u'^: four velocity, and 



illustrate that the results obtained in this scheme ap- 
proximately agree with those in a high-resolution shock- 
capturing scheme based on the Roe-type reconstruction 
for the fluxes [Tsl At each time step, au* is deter- 
mined by solving an algebraic equation derived from the 
normalization u'^u^ = — 1, and then, the primitive vari- 
ables such as p, e, and are updated. 

A uniform atmosphere of small density p « 10'' g/cm'^ 
is added outside neutron stars at i = 0, since the vacuum 
is not allowed in any shock-capturing scheme. However, 
in the HRC scheme, the density can be chosen to be 
much smaller than the previous values |^ . This is a 
benefit in this scheme. The integrated baryon rest-mass 
of the atmosphere is ~ 10~'*Mq in the present simula- 
tion with the largest grid size (see Sec. HI). Hence, the 
effect of the atmosphere for the evolution of binary neu- 
tron stars and its contribution to the disk mass around a 
black hole eventually formed is negligible. (However, for 
the nearly equal-mass merger, the disk mass surround- 
ing a black hole is small as several x10~^Mq; see Sec. 



~dt 



(1) 



For our numerical implementation of the hydrodynamic 
equations, we define a weighted density, a weighted four- 
velocity, and a specific energy defined, respectively, by 
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(2) 
(3) 

(4) 



where h = 1 -I- £ -|- P/ p denotes the specific enthalpy. GR 
hydrodynamic equations are written into the conserva- 
tive form for variables p*, and p*e, and solved using 
a high-resolution central (HRC) scheme 0, . In this 
approach, the transport terms such as di{- ■ ■) are com- 
puted by Kurganov-Tadmor scheme with a third-order 
(piecewise parabolic) spatial interpolation. In |35j . we 



IIV B 31 In such case, we subtract the contribution of the 
atmosphere.) 

As the time slicing condition, an approximate maximal 
slice (AMS) condition K k, Q has been adopted follow- 
ing previous papers 11, 13, 36] to this time. In this 
condition, a is determined by approximately solving an 
elliptic-type equation. This condition is also adopted in 
the case of hypermassive neutron star formation in this 
paper. On the other hand, for the formation of a black 
hole in which an excision evolution is necessary to follow 
its growth, the AMS condition is not advantageous since 
it is not easy to find an appropriate boundary condition 
for a at the excision surface. Thus, for the black hole 
formation, we adopt a dynamical time slicing in which 
the lapse function is determined from 



dtOi = —aK. 



(5) 



In this condition, we need to determine the initial value 
of a. Here, we initially impose the maximal slicing con- 
dition for a. A few simulations were performed in both 
slicing conditions until formation of a black hole, and we 
confirmed that the results depend weakly on the slicing 
conditions, indicating that Eq. jSJ gives a similar slicing 
to the maximal slicing. 

As the spatial gauge condition, we adopt the following 
hyperbolic gauge condition as in [TtIIs^ : 



dtP" =^''\Fi + AtdtF,). 



(6) 



Here, At denotes a time step of the simulation. Successful 
numerical results for the merger of binary neutron stars 
in this gauge condition are presented in [l7l l2]] | . 

In the presence of a black hole, the location is deter- 
mined using an apparent horizon finder for which our 
method is described in The location of the appar- 
ent horizons is used for determining an excision surface 
(see the next subsection). 
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Following previous works, we adopt binary neutron 
stars in close quasiequilibrium circular orbits as the ini- 
tial condition. In computing the quasiequilibrium state, 
we use the so-called conformally flat formalism for the 
Einstein equation (39l |. Solutions in this formalism sat- 
isfy the constraint equations in general relativity, and 
hence, it can be used for the initial condition. The ir- 
rotational velocity field is assumed since it is considered 
to be a good approximation for coalescing binary neu- 
tron stars in nature 4d^. The coupled equations of the 
field and hydrostatic equations ^fj are solved by a spec- 
tral method developed by Bonazzola, Gourgoulhon, and 
Marck Detailed numerical calculations have been 

done by Taniruchi and part of the numerical results are 
presented in |4,1 |. 



B. Black hole excision 

Whenever an apparent horizon is found during the sim- 
ulation, subsequent evolution is followed using an ex- 
cision technique [2^. We adopt the so-called simple 
excision method originally proposed by Alcubierre and 
Briigmann ^2^. In this method, one first determines a 
two-surface inside an apparent horizon which is used for 
the inner boundary in the numerical simulation. And 
then, at the boundary, we impose the so-called copying 
boundary condition in which the time derivative of geo- 
metric variables "fij, (j), Aij, K, and Fi is assumed to be 
spatially constant at the boundary. 

In our method, we choose a cube as the excision bound- 
ary of which the length of each edge is 2£ — 2nAx where 
Aa; denotes the grid spacing, n is a positive integer de- 
termined from the equation 

n = integer[min[rAH(6', ip)]/{V2Ax)] - 1 (7) 

where rAH(^, </3) denotes the coordinate radii of the ap- 
parent horizon. At the surface of the cube, we simply 
copy the time derivative of each variable at the next cell 
along the normal line of each surface. At the corners, we 
copy the time derivative at the next cell along the line 
connecting the corner and origin, and at the edges, at the 
next cell along a perpendicular line connecting the edge 
and another edge in the diagonal direction. 

Since the computational resource is restricted, the radii 
of the apparent horizon are covered only by 6-8 grid 
zones. As a result, the value of n is small; typically n = 2. 
With such a small value, it is difficult to maintain the ac- 
curacy for a long time. However, it is still possible to con- 
tinue the simulation for ~ 0.5-0.7 ms after formation of 
a black hole, which is sufficiently long for approximately 
determining the disk mass around the black hole, and for 
computing ring-down gravitational waveforms associated 
with a black hole quasinormal mode. 



C. Extracting gravitational waves 

Gravitational waves are computed in terms of the 
gauge-invariant Moncrief variables in a flat spacetime [ii * 
as in our series of papers (e.g., IT^^^^|3a. lisf'). The 



detailed equations are describe in jl7l |45j to which the 
reader may refer. In this method, we split the metric in 
the wave zone into the flat background and linear per- 
turbation. Then, the linear part is decomposed using 
the tensor spherical harmonics and gauge-invariant vari- 
ables are constructed for each mode of eigen values {I, m). 
The gauge-invariant variables oil>2 can be regarded as 
gravitational waves in the wave zone, and hence, we fo- 
cus on such modes. In the merger of binary neutron stars 
of mass ratio larger than ^ 0.7, we have found that the 
even-parity mode of (Z, |m|) — (2, 2) is much larger than 
other modes. Thus, in the following, we pay attention 
only to this mode. 

Using the gauge-invariant variables, the luminosity and 
the angular momentum flux of gravitational waves can be 
defined by 



l.m 



dJ 



(8) 
(9) 



Lm 



where i?;™ is the gauge-invariant variable. In this paper, 
we focus only on the even-parity mode with I = 2 for 
Rim- The total radiated energy and angular momentum 
are obtained by the time integration of dE/dt and dJ/dt. 

To search for the characteristic frequencies of gravita- 
tional waves, the Fourier spectra are computed by 



Rlmif) 



Rlm{t)dt, 



(10) 



where / denotes a frequency of gravitational waves. Us- 
ing the Fourier spectrum, the energy power spectrum is 
defined by 



§ = r' E (/>«)' (11) 

i>2,r?j.>0 



where for m 7^ 0, we define 

Rlm{f)=^\Rlm{fW + \Rl-rn{.fW (m>0), (12) 

and use \Rim(,-f)\ = \Rim{f)\ for deriving Eq. 

Since we focus only on the I = 2 even-parity mode, the 
gravitational waveforms are written as 



r 



+ 



— {i?22+(l + cos2 9) cos(2^) 



i?22-(l + cos2 e) sm{2ip)} 

2 , 
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R2or sin 



(13) 
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, 2 

/ix = - 
r 



5 
64^ 



— i?22+ COS 61 s\tl{2lp) 
+ i?22- cos6'cos(2(^) . 
In Eqs. p3|) and (|14ll . the variables are defined by 

R22± 



R22 ± R2 



V2 



In the foUowing, we present 
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-R22-- 



(14) 



(15) 



(16) 



These have the unit of length and provide the ampli- 
tude of a given mode measured by an observer located in 
the most optimistic direction. The amplitude of gravita- 
tional waves, /igw, observed at a distance of r along the 
optimistic direction {9 = 0) is written as 



10 



-22 



i?^ -|- i?x 
0.31 km 



100 Mpc 



(17) 



D. 



Definitions of quantities and metliods for 
calibration 



In numerical simulations, we refer to the total baryon 
rest-mass, the ADM mass, and the angular momentum 
of the system, which are given by 



M 



— I pt:d^x = const. 



(18) 



= — ® diipdSi 

^TT J, — >oo 



(fx, (19) 



1 

8^ 



=60 



d^x, (20) 



where dSj = r'^djrd{cos9)dip, ip^ = —y{dxy + x{dyy , 
Ph = pau^e, Ji — piii, and R/} denotes the Ricci scalar 
with respect to jij . To derive the expressions for M and 
J in the form of volume integral, the Gauss law is used. 

The notations M*i and M*2 are used to denote the 
baryon rest-mass of the primary and secondary neutron 
stars, respectively. In terms of them, the baryon rest- 
mass ratio is defined by Qm = M*2/-M*i(< 1). 

In numerical simulation, M and J are computed us- 
ing the volume integral shown in Eqs. H19|) and (|20|l . 



Since the computational domain is finite, they are not 
constant and decrease after gravitational waves propa- 
gate away from the computational domain. Therefore, 
for t > 0, M and J are not equal to the ADM mass 
and the total angular momentum defined at spatial in- 
finity, but quasilocal quantities. However, in this paper, 
we refer to them simply as the ADM mass and the total 
angular momentum. 

The decrease rates of M and J should be equal to the 
emission rates of the energy and the angular momentum 
by gravitational radiation according to the conservation 
law. Denoting the radiated energy and angular momen- 
tum from the beginning of the simulation to the time t 
as AE{t) and AJ{t), the conservation relations are 



M{t) + AE{t) 
J{t) + AJ{t) = 



= Mo, 
Jo, 



(21) 
(22) 



where Mq and Jq are the initial values of M and J. We 
check that these conservation laws approximately hold 
during the simulation. 

During merger of binary neutron stars (from the last 
one orbit to a relaxed state formed after the merger), 
the angular momentum is dissipated by 15-30% (cf. Sec. 
IVj) . Obviously, the dissipation effect plays a crucial role 
in determining the final outcome. Therefore, checking 
that Eq. H22|l holds in a simulation is one of the most 
important procedures to confirm that the numerical re- 
sults are reliable. 

In addition to checking the conservation of the mass 
and the angular momentum, we monitor the violation 
of the Hamiltonian constraint in the same manner as in 
Ji), 17, 21] in the absence of black hole. In its presence, 
the violation is defined for a region outside the excision 
surface. As demonstrated in 0,1^, the typical mag- 
nitude of the violation is of order 1% throughout the 
simulation in the absence of black hole. In its presence, 
the violation is amplified. However, the typical magni- 
tude is still ~ 10% in the early phase of the excision 
run. The magnitude gradually increases until the crash 
of computation. 



E. Equations of state 

Following |23|, we adopt a hybrid EOS for modeling 
neutron stars' EOS; namely, we write the pressure and 
the specific internal energy in the form 



P = Pcold+Pth, 
£ = Ecold + Eth- 



(23) 
(24) 



Here Pcoid and Ecoid are the cold (zero-tcmpcrature) 
parts, and are written as a function of p. For them, we as- 
sign realistic EOSs for zero-temperature nuclear matter. 
In this paper, we adopt the APR 26] and SLy (Skyrme- 
Lyon) EOSs [S^. These are tabulated as a function of 
the baryon rest-mass density for a wide density range 
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TABLE I: The values of pi in units of c = G = Mq = 1. 
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FIG. 1: Pressure and specific internal energy as a function 
of baryon rest-mass density p for the APR EOS. The solid 
and dotted curves denote the results by fitting formulae and 
numerical data tabulated, respectively. 

from ~ 10 g/cax^ to ^ 10^^ g/cnx^. To simplify numer- 
ical implementation for the simulation, we make fitting 
formulae from the EOS tables as in 'w]. 

In our approach, we first make a fitting formula for 
£coid ill the form 

ecoid(p) = [(1 +Pip^^ +P3/^)(i +P5P^^r - 1] 

X /(-PsP + Pio) 

+ J5l2P^'' /(PSP - Pw)f{-P9P + Pll) 

+ Pl4/"V(P9P-Pll), (25) 

where f{x) = l/(e^ -I- 1), and the constant coefficients 
Pi (i =1-15) are listed in Table I. In making the formulae, 
we focus only on the density for p > 10^° g/cm^, since the 
matter of lower density does not play an important role 
in the merger. Then, the pressure is computed from the 
thermodynamic relation in the zero-temperature limit 

-Pcoid = P — • (26) 
dp 

With this approach, the accuracy of the fitting for the 
pressure is not as good as that in but the first law 



of the thermodynamics is completely satisfied in contrast 
to the work of (46j . 

In Fig. 1, we compare Pcoid and £coid calculated by the 
fitting formulae (solid curves) with the numerical data 
tabulated (dotted curves) for the APR EOS 47). The 
same figures for the SLy EOS are shown in 21]. It is 
found that the fitting formulae agree approximately with 
the tabulated data sets. The relative error between two 
results is within ~ 10% for p > 10^" g/cm^ and less than 
2% for supranuclear density with p ^ 2 x 10^'^ g/cm'^. 

In Fig. [3 we show the relations among the ADM mass 
M, the total baryon rest-mass M*, the central density pc, 
and the circumferential radius R for cold and spherical 
neutron stars in the APR and SLy EOSs. The maxi- 
mum ADM mass (baryon rest-mass) for the APR and 
SLy EOSs is about 2.18 and 2.04Mo, respectively. It is 
worthy to note that for the APR EOS, the radius is in 
a narrow range (11.2-11.4 km) for the ADM mass from 
- 0.5Mq to ~ 1.7Mo, while for the SLy EOS, it de- 
creases with increasing the ADM mass. This results in 
the difference on disk formation in the merger of unequal- 
mass neutron stars. 

Pth and £th in Eqs. (|23|l and (|24|l are the thermal 
(finite-temperature) parts. During the simulation, p and 
e are computed from hydrodynamic variables and e. 
Thus, £th is determined by e — Ecoid- The thermal part of 
the pressure Pth is related to the specific thermal energy 

Sth = £ ~ Ecold as 

Pth = (Fth - l)peth, (27) 

where Tth is an adiabatic constant for which we set Fth = 
2 taking into account the fact that the EOSs for high- 
density nuclear matter are stiff. Since Fth « 5/3 for the 
ideal non-relativistic Fermi gas i* is reasonable to 
consider that it is much larger than 5/3 for the nuclear 
matter. We note that in |23, we also chose Fth = 1-3 
and 1.65 and found that the numerical results depend 
only weakly on its value. 

III. INITIAL CONDITION AND SETTING FOR 
SIMULATION 

In Table II, we summarize several quantities that char- 
acterize irrotational binary neutron stars in quasiequi- 
librium circular orbits used as initial conditions for the 
present simulations. Since the lifetime of binary neutron 
stars from the birth to the merger is longer than ^10^ 
yrs for the observed systems '49^ , the temperature of each 
neutron star will be very low (<^ 10^ K) 50] at the onset 
of merger; i.e., the thermal energy per nucleon is much 
smaller than the Fermi energy of neutrons. Hence, cold 
nuclear EOSs are employed in giving the initial condition. 

We choose binaries of an orbital separation which is 
slightly larger than that for an innermost orbit. Here, 
the innermost orbit is defined as a close orbit for which 
Lagrange points appear at the inner edge of at least one 
of two neutron stars [43. Isif . If the orbital separation 



7 





2.5 


7 (a) 


/ 


^ _ 1 ^ 




2 






- 




1.5 


- 




- 




1 


: 




: 
- 




0.5 




^PR 



















Pc (g/cm3) 





1.5 



0.5 



III 

: APR 


1 1 1 1 1 _ 




(b)- 


y 




1 1 1 1 1 





10 12 
R (km) 



14 



FIG. 2: (a) ADM mass (solid curves) and total baryon rest-mass (dashed curves) as a function of central baryon rest-mass 
density pc and (b) relation between the circumferential radius and the ADM mass for cold and spherical neutron stars in 
equilibrium. "APR" and "SLy" denote the sequences for the APR and SLy EOSs, respectively. 



Model 


M^{Mq) 


Pmax (lO'^g/cm^*) 


Qm 


M,{Mq) 


Mo{Mq) 


go 


Po (ms) 


Co 


Q. 


APR1313 


1.30, 


1.30 


8.62, 8.62 


1.00 


2.858 


2.568 


0.918 


2.064 


0.114 


1.075 


APR1214 


1.20, 


1.40 


8.28, 9.10 


0.842 


2.861 


2.569 


0.920 


2.158 


0.111 


1.076 


APR135135 


1.35, 


1.35 


8.85, 8.85 


1.00 


2.981 


2.665 


0.906 


1.992 


0.120 


1.125 


APR1414 


1.40, 


1.40 


9.09, 9.09 


1.00 


3.106 


2.762 


0.896 


1.923 


0.125 


1.173 


APR1515 


1.50, 


1.50 


9.59, 9.56 


1.00 


3.360 


2.957 


0.879 


1.838 


0.135 


1.269 


APR145155 


1.45, 


1.55 


9.34, 9.86 


0.927 


3.360 


2.959 


0.886 


1.969 


0.129 


1.269 


APR1416 


1.40, 


1.60 


9.09, 10.14 


0.862 


3.363 


2.960 


0.892 


1.969 


0.129 


1.270 


APR135165 


1.35, 


1.65 


8.85, 10.43 


0.800 


3.366 


2.960 


0.888 


1.968 


0.129 


1.271 


APR1317 


1.30, 


1.70 


8.62, 10.74 


0.743 


3.370 


2.960 


0.883 


2.057 


0.126 


1.272 


APR125175 


1.25, 


1.75 


8.40, 11.09 


0.690 


3.377 


2.962 


0.885 


2.145 


0.122 


1.275 


APR1218 


1.20, 


1.80 


8.17, 11.44 


0.639 


3.378 


2.957 


0.861 


2.189 


0.120 


1.275 


SLyl313 


1.30, 


1.30 


8.57, 8.57 


1.00 


2.847 


2.568 


0.922 


2.110 


0.112 


1.175 


SLyl414 


1.40, 


1.40 


9.16, 9.16 


1.00 


3.093 


2.763 


0.902 


2.012 


0.122 


1.277 


SLyl35145 


1.35, 


1.45 


8.85, 9.48 


0.923 


3.094 


2.763 


0.901 


2.013 


0.122 


1.277 


SLyl315 


1.30, 


1.50 


8.56, 9.80 


0.851 


3.096 


2.764 


0.904 


2.104 


0.118 


1.278 


SLyl25155 


1.25, 


1.55 


8.42, 10.15 


0.786 


3.099 


2.765 


0.904 


2.150 


0.117 


1.278 


SLyl216 


1.20, 


1.60 


8.02, 10.54 


0.726 


3.103 


2.765 


0.903 


2.242 


0.113 


1.279 



TABLE II: List of several quantities for initial data (binary neutron stars in quasicircular orbits). The ADM mass of each 
star when they are in isolation A/oo, the maximum density for each star, the baryon rest-mass ratio Qm = M,2/Mti, the 
total baryon rest-mass, the total ADM mass Mo, non-dimensional spin parameter go ~ Jo/M'^, orbital period Po, the orbital 
compactness [Co = (Afof2)^^^], and the ratio of the total baryon rest-mass to the maximum allowed mass for a spherical and 
cold neutron star (Q* = Af*/Aft^niax)- 



becomes smaller than that of the innermost orbit, mass 
transfer sets in and a dumbbell-like structure is formed. 
Until the innermost orbit is reached, the circular orbit 
is stable, and hence, the innermost stable circular orbit 
(ISCO) does not exist outside the innermost orbit for the 
present cases. However, the ISCO seems to be close to 
the innermost orbit since the decrease rates of the energy 
and the angular momentum as a function of the orbital 
separation along the quasiequilibrium sequences are close 



to zero near the innermost orbit. 

The ADM mass of each neutron star, when it is in 
isolation (i.e., when the orbital separation is infinity), is 
denoted by Moo, and chosen in the range between 1.2Mq 
and I.SMq. The labels APR and SLy denote the binary 
models constructed in the APR and SLy EOSs, respec- 
tively. We select the models for a wide range of rest-mass 
ratio, 0.6 <^ Qm < 1, to find the dependence of the disk 
mass around a black hole on Qm for the high-mass case 
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Model 


Grid size 


L (km) 


Ao (km) 


/"merger (kHz) 


Amcrgcr (km 


) Product 


Disk mass (M0) 


APR1313 


(665, 


665, 


333) 


139 


309 


3.20 


94 


NS 


— 




APR1313b 


(377, 


377, 


189) 


78.7 


309 


3.18 


94 


NS 


— 




APR1214 


(377, 


377, 


189) 


74.5 


323 


3.23 


93 


NS 


— 




APR135135 


(377, 


377, 


189) 


77.8 


299 


3.35 


89 


NS 


— 




APR1414 


(665, 


665, 


333) 


135 


288 


3.79 


79 


NS 


— 




APR1515 


(633, 


633, 


317) 


125 


276 


6.5 


46 


BH 


4 X 10" 


4 


APR1515b 


(377, 


377, 


189) 


74.5 


276 






BH 


2 X 10" 


4 


APR145155 


(633, 


633, 


317) 


125 


295 


6.5 


46 


BH 


5 X 10" 


4 


APR1416 


(633, 


633, 


317) 


125 


295 


6.5 


46 


BH 


1.0 X 10^ 


-3 


APR1416b 


(377, 


377, 


189) 


74.5 


295 





_ 


BH 


7 X 10" 


1 


APR135165 


(633, 


633, 


317) 


125 


295 


6.5 


46 


BH 


2.7 X 10" 


-3 


APR135165b 


(377, 


377, 


189) 


74.5 


295 





_ 


BH 


1.9 X 10" 


-3 


APR1317 


(633, 


633, 


317) 


125 


308 


6.5 


46 


BH 


6.9 X 10" 


-3 


APR1317b 


(377, 


377, 


189) 


74.5 


308 


— 


— 


BH 


5.0 X 10" 


-3 


APR125175 


(377, 


377, 


189) 


74.5 


276 





— 


BH 


1.2 X 10" 


-2 


APR1218 


(377, 


377, 


189) 


74.5 


276 





— 


BH 


1.5 X 10" 


-2 


SLyl313 


(633, 


633, 


317) 


131 


316 


3.20 


94 


NS 


— 




SLyl414 


(633, 


633, 


317) 


131 


302 


6.7 


45 


BH 


4 X 10" 


3 


SLyl414b 


(377, 


377, 


189) 


77.8 


302 








BH 


3 X 10" 


3 


SLyl35145 


(377, 


377, 


189) 


77.8 


302 








BH 


5 X 10" 


3 


SLyl315 


(377, 


377, 


189) 


77.8 


315 






BH 


1.4 X 10" 


-2 


SLyl25155 


(633, 


633, 


317) 


131 


322 


6.7 


45 


BH 


4.9 X 10" 


-2 


SLyl25155b 


(377, 


377, 


189) 


77.8 


322 






BH 


3.6 X 10" 


-2 


SLyl216 


(377, 


377, 


189) 


77.8 


336 






BH 


5.7 X 10" 


-2 



TABLE HI: List of setting for simulation and summary of the outcome. L, Ao, and /merger denote the location of outer 
boundaries along each axis, the wave length of gravitational waves at t = 0, and the frequency of gravitational waves from 
the formed hypermassive neutron stars or of ring-down gravitational waves of a quasinormal mode of black holes, respectively. 
Amorgcr deuotes the wave length of gravitational waves Amcrgcr = c/Zmorgcr, which is shown to be L for the large grid sizes. 
In the last two columns, the product at t = 10 ms and disk mass for the black hole formation case are shown. "NS" implies 
that a hypermassive neutron star is the outcome at t = 10 ms, and "BH" implies that a black hole is promptly formed. The 
disk mass is evaluated a.t t — Ikr = 0.5 ms where tAH is the time at apparent horizon formation. Note that we do not compute 
gravitational waves for the small grid size (377, 377, 189) in the case of prompt black hole formation. 



with M ~ 2.76-2.96M0. 

The simulations are performed using a fixed uniform 
grid in the reflection symmetry with respect to the equa- 
torial plane (the orbital plane). In this paper, the typical 
grid size is (633, 633, 317) or (665, 665, 333) for (x, y, z). 
The grid covers the region —L < x < L, —L < y < L, 
and < z < L where L is a constant. The grid spac- 
ing Aa; is ~ 0.4 km with which the major diameter of 
each star is covered with about 45 grid points initially. 
We have shown that with this grid spacing, a convergent 
numerical result is obtained in the merger simulations 

m 

Accuracy in the computation of gravitational wave- 
forms (in particular amplitude) and the radiation reac- 
tion depends on the location of the outer boundaries if 
the wavelength, A, is larger than L For L ^ 0.4A, 

the amplitude and the radiation reaction of g ravitational 
waves are significantly overestimated [T7ll52j| Due to the 
restriction of the computational power, it is difficult to 



take a huge grid size in which L is much larger than A. 
As a consequence of the present restricted computational 
resources, L has to be chosen as ~ 0.4-0.45Ao where Aq 
denotes the value of A at i = 0. Hence, the error asso- 
ciated with the small value of L is inevitable; the am- 
plitude and radiation reaction of gravitational waves are 
overestimated in the early phase of the simulation. The 
overestimation of the radiation reaction leads to slight 
spurious shortening of the late inspiraling phase. 

However, the typical wavelength of gravitational waves 
just before the merger quickly becomes shorter due to 
a short radiation-reaction time scale, and hence, the 
accuracy of the wave extraction is improved with the 
evolution of the system. In particular, the wavelength 
of quasiperiodic gravitational waves emitted from the 
formed hypermassive neutron star and of ring-down grav- 
itational waves associated with a quasinormal mode of 
the formed black hole (denoted by Amerger) is much 
shorter than Aq and satisfies the condition Amcrgcr < L 
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for the grid size (633, 633, 317) (see Table III). There- 
fore, the waveforms in the merger stage are computed 
accurately (within ~ 10% error) as confirmed in p^Tj . 

For some models, the simulations are performed for a 
smaller grid size (377, 377, 189) but with the same grid 
spacing Aa;. In this case, gravitational waveforms cannot 
be accurately computed for the inspiraling phase because 
of the small value of L/Xq. However, if we pay atten- 
tion only to the outcome after the merger, such small 
value for L may be allowed since the properties of the 
outcome depend weakly on L. This fact is indeed con- 
firmed in the simulations for several models (cf. Table 
III; about the reason for the systematic underestimation 
of the disk mass, see IIV CTTl . Thus, we perform simu- 
lations for models APR135135, APR1214, APR125175, 
APR1218, SLyl35145, SLyl315, and SLyl216 only with 
the small grid size. In this case, the focus is only on the 
outcome, not on computing gravitational waveforms. 

With the (633, 633, 317) grid size, about 240 GBytes 
computational memory is required. For the case of the 
hypermassive neutron star formation, the simulations are 
performed for about 30,000 time steps (until i ~ 10 ms) 
and then stopped to save the computational time. The 
computational time for one model in such a simulation 
is about 180 CPU hours using 32 processors on FACOM 
VPP5000 in the data processing center of National As- 
tronomical Observatory of Japan (NAOJ). For the case 
of the black hole formation, the simulations are contin- 
ued until mass accretion rate from the surrounding disk 
to a formed black hole relaxes approximately to a con- 
stant. In this case, the computational time is about 60 
CPU hours for about 12,000 time steps. 



IV. NUMERICAL RESULTS: MERGER 
PROCESS 

A. Summary of the outcome 

For all the cases, the binary orbit is stable at t — 0, 
but after slight decrease of the orbital separation due 
to gravitational radiation reaction, merger sets in. In 
the present simulations, the merger starts in one orbit 
{in t ^ 2 ms) irrespective of the models. If the to- 
tal ADM mass of the system is high enough, a black 
hole is formed within ^ 1 ms after two stars come into 
contact. On the other hand, for models with the ADM 
mass smaller than a threshold mass Mthr, a hypermas- 
sive neutron star is formed and survives for more than 
10 ms. However, it will collapse to a black hole even- 
tually due to gravitational radiation reaction or angu- 
lar momentum transport by other effects (see discussion 
in Sec. IIV B 31) . In the case of black hole formation 
(models APR1515, APR1416, APR135165, APR1317, 
APR125175, APR1218, SLyl414, SLyl35145, SLyl315, 
SLyl25155, and SLyl216), the evolution of the black 
holes is followed using a black hole excision technique 
I29I until the mass accretion rate from the disk to the 



formed black hole becomes approximately constant. 
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FIG. 3: Evolution of the central value of the lapse func- 
tion ttc and the maximum values of the rest-mass density 
Pmax for models APR1313 (solid curves), APR135135 (long- 
dashed curves), APR1414 (dashed curves), and APR1515 
(dot-dashed curves). The dotted horizontal lines denote the 
central values of the lapse and rest-mass density for the 
marginally-stable, spherical neutron star in equilibrium with 
the cold APR EOS. 

In Fig. 3, we show the evolution of the central 
value of the lapse function ac and the maximum of 
the baryon rest-mass density pmax for models APR1313, 
APR135135, APR1414, and APR1515. For model 
APR1515, Uc collapses to zero and Pmax quickly increases 
at i ^ 1.5 ms, implying that a black hole is formed 
promptly soon after the onset of the merger (at t « 2.04 
ms). For other cases, ac and Pmax settle down to relaxed 
values, implying that a quasistationary neutron star is 
formed. Since the mass is larger than the maximum al- 
lowed limit of rigidly rotating neutron stars, these neu- 
tron stars are hypermassive |23l |. 

For model APR1414, ac and Pmax oscillate with an 
amplitude larger than that for models APR1313 and 
APR135135. This behavior results from the fact that 
the self-gravity is large enough for the merged object to 
deeply shrink surmounting the centrifugal force. This in- 
dicates that the total ADM mass of this model (M « 
2.76Mq) is slightly smaller than the threshold value for 
the prompt black hole formation. Since a black hole is 
formed for M « 2.96M0, the threshold mass for black 
hole formation is Mthr ~ 2.8-2.9M0. In a previous 
paper |2l|, we adopted the SLy and FPS (Friedman- 
Pandharipande-Skyrme) EOSs and found that Afthr ~ 
2.7 Mq and 2.5Mq, respectively. The high value of Afthr 
for the APR EOS is reasonable since it is stiffer than 
the SLy and FPS EOSs, and the maximum allowed mass 
for spherical, cold neutron star Msph is also larger. The 
present result combined with the previous one |2l| sug- 
gests an empirical relation Mthr/A^sph = 1.30-1.35 for 
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these stiff nuclear EOSs. 

We note that this relation is highly different from that 
for the r-law EOS with T = 2 |[17(| for which Mthr/Mgph « 
1.7. As discussed in the compactness of each neutron 
star in the stiff nuclear EOSs is larger than that with the 
r = 2 EOS for a given mass. Accordingly, for a given 
total mass, the binary system at the onset of the merger is 
more compact. This implies that the angular momentum 
is dissipated more before the merger sets in with the stiff 
nuclear EOSs. The dissipation of the angular momentum 
by ~ O.lJo before the merger helps the prompt black hole 
formation, since the hypermassive neutron star formation 
requires a substantial centrifugal force to sustain the self- 
gravity 

In recent papers 53] , ratio of the maximum mass of hy- 
permassive neutron stars in equilibrium to Mgph has been 
investigated in detail. These works coincidently show a 
similar dependence of the ratio Mthr/A/gph on the EOSs. 
Hence, the small value of Mthr/A^sph may be partly due 
to the absence of high-mass differentially rotating neu- 
tron stars in equilibrium. 

In Figs. 0HH1 display the snapshots of the density 
contour curves and the velocity vectors in the equato- 
rial plane at selected time slices for models APR1313, 
APR1414, APR1515, APR1416, and APR1317, respec- 
tively. In the first two hypermassive neutron 
star is formed, while a black hole is a prompt outcome in 
other three cases. For models APR1515, APR1416, and 
APR1317, the ADM mass of the system is approximately 
identical while the mass ratio is different. The struc- 
ture and density of disk surrounding the formed black 
hole depend significantly on the mass ratio as found from 
Figs. EllHl Figure ini displays the density contour curves 
and the velocity vectors in the y = plane for models 
APR1313, APR1414, APR1515, APR1416, APR135165, 
and APR1317 at a late time when the mass accretion rate 
relaxes approximately to a constant. This shows that (i) 
the hypermassive neutron stars have a highly flattened 
structure and (ii) the disk surrounding the black hole is 
geometrically thin for the mass ratio close to unity but 
can be thick for the smaller mass ratios with Qm ^ 0.8. 

Figure^lshows the angular velocity along the x and y 
axes (a) of relaxed hypermassive neutron stars for models 
APR1313 and APR1414 and (b) of the disk surrounding 
the black hole for model APR135165. This illustrates 
that the hypermassive neutron stars, in particular in their 
outer region, are differentially rotating and the central 
part is rapidly rotating with the rotational period shorter 
than 1 ms. The disk surrounding the black hole has ap- 
proximately a Kepler orbit outside the ISCO (for which 
the estimated coordinate radius is ^ 10 km). 

In the following two subsections, details about the for- 
mation process of hypermassive neutron stars and black 
holes are discussed separately. Implication of the results 
to formation of a central engine of SGRBs is also dis- 
cussed in the subsequent section HVDI 



B. Formation of hypermassive neutron star 

1. Models APR1313 and APR135135 

In the formation of the hypermassive neutron stars 
from equal-mass neutron stars with relatively low total 
mass ~ 2.6-2. 7Mq (models APR1313 and APR135135), 
a double core structure is first formed (see the snapshot 
aXt = 3.218 ms of Fig. 0)), and then, it relaxes to a highly 
nonaxisymmetric ellipsoid [see the snapshots for i > 4 ms 
of Fig. 2|and|2a)]. The contour plots, drawn for a high- 
density region with the rest-mass density larger than the 
nuclear density ~ 2 x 10"'^^ g/cvc?, show that the axial 
ratio of the ellipsoid measured in the equatorial plane is 
< 0.8 for i ^ 4 ms. Figure|^a) also shows that the axial 
length along the z-axis for p > 2 x 10^* g/cm^ is about 
7 km, which is ~ 70% of the major axis. Namely, a ro- 
tating ellipsoid of a large ellipticity is the outcome. This 
result is essentially the same as that found in '2l'| for 
models SLyl313 and SLyl25135. This large ellipticity is 
achieved due to the rapid rotation and the high stiffness 
of the chosen EOS. 

The rapid rotation is found from Fig. Iior al: it shows 
that the rotational period of the central region is shorter 
than 1 ms at i ~ 10 ms, at which the hypermassive 
neutron star relaxes approximately to a stationary state. 
This indicates that the rotational centrifugal force plays 
an important role for sustaining the large self-gravity of 
the hypermassive neutron stars [2^ ■ 

For the unequal-mass case with Qm ~ 0.85 (model 
APR1214), the smaller-mass neutron star is tidally elon- 
gated at the merger. As a result, the double-core struc- 
ture is not formed in contrast to model APR1313. How- 
ever, the formed hypermassive neutron star relaxes to an 
ellipsoid in the similar way to model APR1313. Subse- 
quent evolution proceeds in essentially the same manner 
to that in the equal-mass case. 

The merger for model APR135135 proceeds in qualita- 
tively the same way as that for model APR1313. Quan- 
titatively, the results such as the ellipticity and the ro- 
tational period of the hypermassive neutron star are 
slightly different reflecting the difference in the mass; el- 
lipticity is slightly smaller and the rotational period is 
slightly shorter because of its larger compactness. 

Because of the nonaxisymmetric structure, the hy- 
permassive neutron stars found for models APR1313, 
APR1214, and APR135135 emit quasiperiodic gravita- 
tional waves and the angular momentum is dissipated 
substantially. However, the dissipation time scale is much 
longer than 10 ms, implying that they remain the ellip- 
soidal star for longer than 10 ms (cf. Sec. 0. These ellip- 
soids are also differentially rotating (cf. Fig. I10|) . Thus, 
they are subject to angular momentum transport by mag- 
netic effects |23|5^|5^0. However, note that the time 
scale of the angular momentum transport is shorter than 
~ 100 ms only when very strong magnetic fields greater 
than ^ 10^^ G are present p^ . 

Assume that the dissipation time scale by gravitational 



11 




20 10 10 20 
X (km) 



10 10 20 
X (km) 



FIG. 4: Snapshots of the density contour curves for p in the equatorial plane for model APR1313. The solid contour curves are 
drawn for p = 2 x lO" x i g/cm^ {i = 1, 2, 3, • • •) and for 1 x IQ^''-''-^' g/cm^ (i = 1 ~ 6). The (blue) thick dotted and solid 
curves denote 1 x 10^* g/cm^ and 1 x lO^'^ g/cm^, respectively. The number in the upper left-hand side denotes the elapsed 
time from the beginning of the simulation in units of ms. Vectors indicate the local velocity field {v^,v^), and the scale is 
shown in the upper right-hand corner. 



waves is shortest among other processes. Then, there 
are two possible fates after a longterm emission of grav- 
itational waves. One fate is that after the angular mo- 
mentum dissipation, the centrifugal force becomes weak 
enough to induce gravitational collapse to a black hole. 
The other is that the angular momentum (i.e., rotational 
kinetic energy) is dissipated to be too small to main- 
tain the nonaxisymmetric structure and, consequently, a 
spheroidal star in a stationary state is formed. (It is well 
known that an elliptical structure can be achieved only 
in the case that the ratio of the rotational kinetic en- 
ergy to the gravitational binding energy is large enough 
[53 ■) If the time scale for the decrease of the angular 
momentum is shorter (longer) than that for the decrease 
of the ellipticity, a black hole (a spheroidal hypermas- 
sive neutron star) is the outcome. In the formation of 
a spheroid, subsequent evolution will be determined by 
other processes such as angular momentum transport by 



the magnetic braking [23| or the magnetorotational in- 
stability (MRI) 55]. In the previous paper in which 
the SLy and FPS EOSs are adopted, we have found that 
a black hole is the outcome (see the results for models 
SLyl35135b and FPS125125b in HH). As the merger 
process depends weakly on the EOSs, the black hole for- 
mation may be also the fate in the APR EOS. However, 
this may not hold as illustrated in the next section for 
high-mass model APR1414. To answer this question, an 
extremely longterm simulation up to t ~ 100 ms is nec- 
essary. This problem is left as the future issue. 



2. Model APRI414 

In the formation of the hypermassive neutron stars 
with high mass 2.8Mq (model APR1414), the evolu- 
tion proceeds in a different manner from that for models 
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FIG. 5: The same as Fig. glbut for model APR1414. 
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APR1313 and APR135135. Since the mass is larger, a 
more compact merged object is formed soon after the 
onset of the merger (second panel of Fig. |3J). Then, it 
bounces back to a state with a large radius (or a small 
density; cf. Fig. 3) and repeats quasiradial oscillations 
with a large amplitude before relaxing to a hypermas- 



sive neutron star in a quasistationary state. At the state 
of the large radius, the angular momentum is efficiently 
transported from the inner region to the outer envelop 
(see discussion below). Oir the other hand, the compact 
state achieved at maximum compressions (cf. second and 
fourth panels of Fig. |SJ) results in a large amount of angu- 
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FIG. 8: The same as Fig. Elbut for model APR1317. 



lar momentum dissipation via gravitational radiation (cf. 
SecEJ- These effects subsequently lead the hypermas- 
sive neutron star to a state with relatively small angular 
momentum. 

The decrease of the angular momentum in a shorter 
time results in a state of small rotational kinetic energy. 
On the other hand, the gravitational binding energy is 
larger than in models APR1313 and APR135135, be- 
cause of its larger compactness. As mentioned in the 
last part of Sec. IIVB II an elliptical structure can be 
achieved only in the case that the ratio of the rotational 
kinetic energy to the gravitational binding energy is large 
enough '58*1 . For the hypermassive neutron star formed in 
model APR1414, this ratio should be smaller than those 
in models APR1313 and APR135135. Consequently, a 
spheroid of small ellipticity is formed [cf. the last three 



panels of Fig. |S1 and IHIb)] . This results in the facts that 
luminosity of gravitational waves is much smaller than 
for model APR1313 for i ^ 5 ms (cf. Sec. |V||, and the 
dissipation time scale of the angular momentum by grav- 
itational radiation is much longer. In model APR1414, 
probably, dissipation by gravitational radiation will not 
induce gravitational collapse. 

As in the hypermassive neutron star for model 
APR1313, the outer region with the coordinate radius 
i? <; 8 km is differentially rotating [see Fig. llUf a')]. Such 
differentially rotating region is subject to the magnetic 
braking and the MRI in the presence of magnetic fields 
[ssl ISa . Is^ . If the instabilities turn on, the magnetic 
fields play an important role for transporting angular 
momentum outward, subtracting the angular momentum 
of the central part. Hence, a plausible scenario for the 
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t = 2.959 ms for APR135165, and (f) at t — 3.038 ms for APR1317. The contour curves and the velocity vectors are drawn in 
the same way as in Fig. El 



hypermassive neutron star of model APR1414 is the col- 
lapse to a black hole by the magnetic effects. 

Finally, we note that in the previous paper ,21.] in 
vifhich the SLy and FPS EOSs are used, we have not found 
the formation of spheroids, but rather, gravitational radi- 
ation triggers the collapse of hypermassive neutron stars 
to a black hole. The likely reason in this difference is 
that the APR EOS is stiffer than others. Namely, the 
hypermassive neutron stars can be more compact with 
this EOS [cf. Fig. ^h)] escaping the collapse. In such 
extremely compact state, the gravitational binding en- 
ergy would be large enough to reduce the ratio of the 
rotational kinetic energy to the binding energy below the 
threshold value of the formation of an ellipsoid. 



3. Angular momentum transport and disk mass 

Due to a torque by the nonaxisymmetric structure of 
the merged object, angular momentum is transported 
from the inner region to the outer one. As a result, 
the rotational angular velocity fl = v"^ near the center 
decreases and its profile is modified. Figure ^1 shows 
the angular velocity Q of the hypermassive neutron stars 



along X and y axes at t ^ 10 ms. The hypermassive 
neutron stars are differentially and rapidly rotating at 
the birth [23|; the rotational period around the central 
region is « 0.2-0.3 ms at the birth. (The rotational pe- 
riod of the double core found in Figs. 01 and [S] which is 
formed soon after the merger, is ^ 0.2-0.3 ms. This is 
found from the frequency of gravitational waves; cf. Sec. 
IVj) . The angular momentum is then redistributed by the 
torque, resulting in a fairly uniform profile of fl near the 
rotational axis with R <, 8 km. 

The significance of the angular momentum transport 
is also found from the evolution of the mass spectrum 
as a function of the specific angular momentum, M^{j). 
Here, the specific angular momentum j is defined by huip 
and the mass spectrum M, (jo) is given by an integrated 
baryon rest-mass of fluid elements with j < jo; 

M^Jo) = f P*<i^x. (28) 

j<io 

In Fig. m we show this mass spectrum at selected 
time slices for models APR1313 and APR1414. This 
illustrates that the fraction of baryon rest-mass with 
j ^ 1.4Mo is absent at t = but increases with time. 
This is due to the angular momentum transport from the 
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central hypermassive neutron star of ellipsoidal shape to 
fluid elements in the outer envelop by the torque associ- 
ated with the nonaxisymmetric structure. However, the 
fraction of the rest-mass with j > 2Mo is still ~ 1% 
of the total mass because such fluid elements are absent 
initially, and also probably because of gravitational ra- 
diation which carries away the angular momentum by 
~ 30% in 10 ms (cf. Sec. 

Assume a hypothetical case in which the central region 
collapses to form a black hole at t w 10 ms. Then, the 
ADM mass and the angular momentum of the black hole 
will be « 0.97Mo and ~ O.TJq, respectively, implying 
that the spin parameter is g ~ 0.7 (cf. Sec. 0). The 
specific angular momentum at the ISCO, iiscm around a 
black hole with q = 0.7 is « 2.58M [sg as plotted by 
the dotted vertical line in Fig. ^] The mass of the fluid 



elements with j > 2.5M is « O.OIMq at t = 10 ms for 
model APR1313. Namely, a disk of mass ~ O.OIM© will 
be formed in this hypothesis. The value for the disk mass 
is approximately identical with that presented in [sij in 
which a similar nuclear EOS 32] is adopted. 

We note that the lifetime of the hypermassive neutron 
star is in reality much longer than 10 ms (assuming no 
other strong dissipation process than gravitational radi- 
ation), and hence, the angular momentum will be dissi- 
pated much more by gravitational radiation before for- 
mation of a black hole. Therefore, the value of q will 
be smaller and jisco could be larger (e.g., for q = 0.5, 
jisco ~ 2.90M). This effect will slightly reduce the disk 
mass. On the other hand, the angular momentum will 
be gradually transported in the outer region with time 
due to the torque from the central region of a nonaxisym- 
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metric shape. This effect wiU shghtly increase the disk 
mass. 

For model APR1414, the fraction of the rest-mass with 
j > 2.5M is slightly larger as « O.O2M0 at t R:! 10 ms. 
The reason is that in this case, the merged object quasira- 
dially oscillates with a larger amplitude soon after the on- 
set of the merger (cf. Fig. . During the high-amplitude 
oscillation, matter in the outer region expands to a large 
radius and can gain a torque from the ellipsoidal hyper- 
massive neutron star for a longer time. Actually, the 
fraction of the rest-mass with j > 2.5M steeply increases 
for 4 ms ^ t <, 6 ms during which the hypermassive 
neutron star repeats quasiradial oscillations of the high 
amplitude. The high-amplitude oscillation results from 
the fact that the ADM mass of the system is close to 
Mthr- Hence, for the system of M ^ Mthr, the angular 
momentum transport efficiently works. 

The outcome for model APR1414 at < ~ 10 ms is a hy- 
permassive neutron star of nearly spheroidal shape. This 
implies that the angular momentum will not be signif- 
icantly dissipated any longer and also outward angular 
momentum transport will not work efhciently for t > 10 
ms. Therefore, in this case, the disk of mass at least 
« 0.02Mq will be formed if the hypermassive neutron 
star collapses to a black hole by other mechanisms. 

For model APR1214, for which the ADM mass is ap- 
proximately identical with that for model APR1313, the 
less massive star is tidally elongated by the primary at the 
merger and subsequently forms an accretion disk around 
a formed hypermassive neutron star (evolution of the 
contour curves is similar to those of model APR1416 for 
t <, 2 ms; cf. Fig. [7|). Because of the nonaxisymmet- 
ric density profile, the angular momentum is transported 
outward more efficiently than in model APR1313. How- 
ever, the baryon rest-mass with j > 2.5M for model 
APR1214 at t = 10 ms is « O.O25M0 which is only 
slightly larger than that for model APR1313, and hence, 
the degree of the increase is not as drastic as that re- 
ported in jSlj] for Qm ^ 0.85. The likely reasons for our 
small mass are (i) the radius of neutron stars depends 
very weakly on each ADM mass in the APR EOS (cf. 
Fig. |2Jl and (ii) gravitational radiation carries a large 
amount of angular momentum (cf. Sec. which is 
not taken into account in the work of [sj- Because of 
the weak dependence of the stellar radius on the mass, 
the tidal deformation becomes important only for close 
orbits. At such small orbital radius, the angular momen- 
tum is already dissipated by gravitational radiation to 
be small, and hence, the angular momentum transport 
by a torque associated with the nonaxisymmetric struc- 
ture does not help sufficiently increasing the mass with a 
large value of j . 

Formation of hypermassive neutron stars from binaries 
of small values of Qm ^ 0.8 is unlikely since the low- 
mass binary neutron star with such small value of Q]\i 
would be absent; i.e., the mass of neutron stars should 
be larger than ~ 1.2Mq while the total mass has to be 
smaller than Mthr, implying that Qm should be larger 



than 0.8. Therefore, the fraction of the mass with 
j ^ 2.5M will be at most ~ O.OSM© « O.OlAf, for the 
hypermassive neutron star formation case in the APR 
EOS. 

For model SLyl313 in which a hypermassive neutron 
star is also formed, the merger process and subsequent 
evolution of the hypermassive neutron star is quite sim- 
ilar to those for models APR1313 and APR135135 
As a consequence, the rest-mass with j > 2.5M at t = 10 
ms is the similar value ~ 0.02Mo. This suggests that 
the value of the disk mass depends weakly on the EOS. 
However, we note that the disk mass could be larger if 
we employ EOSs with which the neutron star radius is 
larger than that of the APR and SLy EOSs, and hence, 
the gravitational radiation reaction at the merger is less 
important. 

Before closing this section, we note that the discus- 
sion about the disk mass here is based on the assump- 
tion that angular momentum is transported only by a 
torque in the hypermassive neutron stars. In the pres- 
ence of magnetic fields with a sufficient strength, the 
angular momentum could be transported efficiently in 
a short time scale |i56. 57.] ■ and as a result, the hypermas- 
sive neutron stars such as formed for models APR1313, 
APR1414, APR1214, and SLyl313 may collapse ejecting 
a large amount of the mass to form a disk. The system 
eventually formed could be composed of a black hole and 
a massive disk of ^ 0.05Mq as illustrated in the recent 
magnetohydrodynamic simulations [56lf57| . 



C. Formation of black hole 

For binary neutron stars with M w 2.96Mo in the 
APR EOS and M « 2.76Mo in the SLy EOS, a black 
hole is formed soon after the onset of the merger irre- 
spective of the mass ratio Qm- Until two stars come into 
contact, evolution proceeds in the similar manner to that 
in the formation of hypermassive neutron stars. However, 
after the contact, the merged objects quickly collapse to a 
black hole. For the unequal-mass case, a less massive star 
is tidally deformed just before the merger. The degree of 
the deformation as well as lag angle, which is defined to 
be the angle in the equatorial plane between the major 
axis of each star and the axis connecting the centers of 
mass of two stars, are larger for smaller values of Qm- 
Because of the nonzero lag angle, angular momentum 
is transported outward by a torque from the nonaxisym- 
metric merged object, increasing the specific angular mo- 
mentum of the matter located in the outer region. Only 
in the case that the angular momentum transport works 
efficiently, an accretion disk may be formed around the 
central black hole since the specific angular momentum 
before the onset of the merger is too small as illustrated 
in Sec. IIV B 31 fsee also a discussion below). 
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FIG. 12: (a) Evolution of the baryon rest-mass of disks sur- 
rounding black holes and (b) the evolution of irreducible mass 
of black holes for models APR1515 (solid curve), APR145155 
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SLyl25155 (dotted curve), ^ah denotes the time at the first 
formation of an apparent horizon. 



0.1 



0.01 



w 
w 



w 0.001 



Q 



0.0001 



r 1 1 1 1 1 




II 




□ 




- 




- 

5i2.76MQ - 














• ^ s 














• \ 






^ -o 








" - f 


1 1 1 1 1 1 1 


APR, M 





0.7 



FIG. 13: Baryon rest-mass of disks around a black hole as 
a function oi Qm = M*2/M,i for the black hole formation 
cases. The disk mass is evaluated at t — tAn = 0.5 ms. The 
open circles and squares denote the results with (633, 633, 
317) grid size in the APR and SLy EOSs, respectively. The 
filled circles and squares denote the results with (377, 377, 
189) grid size in the APR and SLy EOSs. The dashed curves 
denote the fitting formulae (I29I I. The results for the APR 
and SLy EOSs are derived for M f=i 2.96Mq and 2.76Mo, 
respectively. 



In Fig. I12f a). we plot the total baryon rest-mass lo- 
cated outside the apparent horizon for models APR1515, 
APR145155, APR1416, APR135165, APR1317, and 
SLyl25155 as a function of t — tAH where iAH denotes 
the time at the first formation of an apparent horizon. 
In Fig. ^1 the baryon rest-mass at t — ^ah = 0.5 ms as a 
function of Q m is also plotted for all the models in which 
a black hole is formed. These figures show that the final 
state is composed of a central black hole and a disk of a 
small mass. A rapid accretion proceeds in the first ~ 0.2 
ms after formation of apparent horizon. In particular, 
for Qm ^ 0.9, most of the fluid elements are swallowed 
into the black hole in such a short time. With decreasing 
the value of Qm, the fraction of the fluid elements which 
escape falling into the black hole steeply increases. 

We note that for the smaller grid size, the disk mass is 
smaller even for the same model (cf. Table III) . As men- 
tioned in Sec. IIIII the inspiral time is spuriously shorten 
with the smaller grid size since the radiation reaction is 
overestimated due to the situation with L ^ Aq. As 
a result, the time for transporting angular momentum 
outward is shorten, resulting in the decrease of the disk 
mass. 

For Qm ^ 0.9, the baryon rest-mass of the disk be- 
comes smaller than O.OIMq within ~ 0.3 ms after for- 
mation of apparent horizon (« 5 x 1O~^M0 for models 
APR1515 and APR145155, and w 5 x IQ-^Mq for mod- 
els SLyl414 and SLyl35145). The fundamental reason 
is that the specific angular momentum j for all the fluid 



elements before the merger is too small. In the present 
case, the maximum value of j in the initial conditions 
is ~ 1.3M, which is much smaller than the value of j 
at the ISCO of a Kerr black hole of mass M and spin 
parameter 0.7-0.9. Moreover, the specific angular mo- 
mentum decreases due to emission of gravitational waves. 
Hence, a disk can be formed only in the case (i) an efh- 
cient mechanism for angular momentum transport which 
works within ~ 1 ms after the onset of the merger (before 
collapsing to a black hole) is present and (ii) the trans- 
port mechanism works efficiently enough for overcoming 
the dissipation by gravitational radiation in the outer re- 
gion. For Qm 1, angular momentum transport due to 
a torque by the nonaxisymmetric central object does not 
work efficiently, resulting in the small disk mass. 

We note that in the prompt formation of a black hole, 
the time for the angular momentum transport is much 
shorter than in the formation of a hypermassive neutron 
star. This is the main reason that the disk mass is much 
smaller than that in Sec. IIV B 31 for nearly equal-mass 
binaries. 

For Qm ^ 0.9, the disk mass increases with decreas- 
ing the value of Qm because the less massive neutron 
star is tidally elongated at the merger, and hence, efh- 
ciency of the angular momentum transport is enhanced. 
However, the disk mass does not exceed O.OIMq even for 
Qm = 0.75 for the APR EOS. The hkely reason is that 
the stellar radius is small irrespective of the mass of each 
star in this EOS. Due to this property, the merged object 
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is compact at the onset of merger enough to promptly 
produce a black hole without transporting angular mo- 
mentum efficiently. 

The disk mass is much larger in the SLy EOS for a 
given value of Qm than in the APR EOS. The reasons are 
(i) the mass (M « 2.7QMq) is smaller than for the APR 
EOS case (M « 2.96M0) and (ii) with the SLy EOS, 
the stellar radius of relatively small mass M ^ 1.3Mq 
is larger (cf. Fig. With the smaller mass, the grav- 
itational radiation reaction time scale is longer. Con- 
sequently, the less massive star in the binary is tidally 
elongated by a larger degree, slightly postponing the col- 
lapse of the merged object to a black hole and enhancing 
efficiency of the angular momentum transport. 

It is interesting to note that the disk mass decreases 
during the growth of a black hole in the following mecha- 
nism: In the merger of unequal-mass binary, the less mas- 
sive neutron star is always tidally elongated. The inner 
part first collides the companion, but slips through the 
companion's surface to form a small spiral arm around 
a black hole which is subsequently born. Outer part of 
larger angular momentum also forms a large spiral arm. 
The phase difference between these two spiral arms is 
about 180° at their formation (cf. Figs. EHHl). However, 
since the smaller spiral arm rotates faster than the larger 
one, the two spiral arms collide each other to generate 
shocks in one rotational period, which convert a large 
fraction of the kinetic energy to the thermal energy. As 
a result, an amount of matter is swallowed into the black 
hole quickly. The collision of the spiral arms also leads 
to geometrical thickening of the disk [cf. Fig. (SJe) and 
(f)] as well as to the heat-up. 

For M « 2.96M0 in the APR EOS and M w 2.76Mq 
in the SLy EOS, the numerical results of the disk mass 
Md as a function of Qm for 0.75 ^ Qm < 1 is approxi- 
mately fitted by 

Md^Mdo + Mdiil-QMf, (29) 

where Mao = 0.0004Mq, Mai ~ I.UMq, and p = 4 for 
M « 2.96M0 in the APR EOS with (633, 633, 317) grid 
size, and Mao = O.OO3M0, Mdi « 3.33Af0, and p = 3 for 
M « 2.76M0 in the SLy EOS with (377, 377, 189) grid 
size. Since M^o is small, Md is approximately propor- 
tional to (1 — QmY for Qm ^ 0.8, and hence, the value 
of Mdi essentially determines the relation. As mentioned 
previously, the value of Md varies by ~ 30% depending 
on the grid size by which the radiation reaction time scale 
and resulting transport time of the angular momentum 
spuriously change. 

As indicated by this fitting formulae (|29|l . the disk 
mass steeply increases with decreasing the value of Qm 
(cf. Fig. 113(1 , although for Qm ^ 0.9, the disk mass is 
much smaller than 0.01 A/0. For Qm ^ 0.7, the increase 
rate of Md with decreasing value of Qm is not as steep 
as Eq. (^5)) . but it is still positive. 

We note that the disk mass obtained here for the 
prompt black hole formation case gives an approximate 
upper limit for both the APR and SLy EOSs, since the 



ADM mass adopted is close to Mthr- For more massive 
case, the system at the merger becomes more compact 
and has a shorter dynamical time scale. Therefore, out- 
ward transport of angular momentum works for a shorter 
time scale, resulting in smaller disk mass. 

Figure [rUT b) shows the angular velocity of the disk for 
model APR135165 and illustrates that the disk approx- 
imately has the Kepler angular velocity for the radius 
larger than ^ 10 km which is approximately the radius 
of the ISCO around the black hole. The angular velocity 
of the inner part with i? ^ 10 km is not as large as the 
Kepler one, indicating that matter in such region grad- 
ually falls into the central black hole. This is also seen 
from the velocity fields in Fig. injd)-(f). In reality, the 
matter in the accretion disk of the Kepler angular veloc- 
ity will subsequently evolve in a dissipation time scale, 
which is determined either by viscosity or by magnetic 
fields. 

For Qm ^ 0.9, not only the disk mass is small but also 
the disk is geometrically thin. These properties are com- 
pletely unfavored for producing high-Lorentz factor jets 
of GRBs 6]. Therefore, prompt formation of a black hole 
in the merger of binary neutron stars of the mass ratio 
Qm ^ 0.9 is unlikely to be a scenario for producing a 
central engine of SGRBs. On the other hand, for a suffi- 
ciently smaU value of Qa/ [Qm ^ 0.75 for the APR EOS 
and Qm ^ 0.85 for the SLy EOS), the disk mass will be 
larger than ~ O.OlAf0. In addition, the disk is geometri- 
cally thick as indicated in Fig. for model APR1317; 
i.e., a torus is formed. Furthermore, the thermal energy, 
which is generated in the shock heating between spiral 
arms in the disk and is evaluated from e — ecoid PH , is 
high enough (typically ~ 1-2 x 10^^ K) for generating a 
large amount of thermal neutrinos \wL \6U^ . Viscous heat- 
ing will also play a role for subsequently keeping such 
high temperature. Implication of such hot torus with 
mass J; O.OlAf0 will be given in Sec. rrVDl 



2. Properties of black hole 

The area of the apparent horizon Aah is determined 
in the black hole formation cases. From the area Aah, 
the irreducible mass may be approximately defined by 

M..,.^. (30) 

which varies from 0.9Mo to 0.97A/o for all the models in 
the APR EOS as shown in Fig. [Elb). Since most of 
the fluid elements are swallowed into the black hole and, 
also, the energy carried out by gravitational radiation is 
w O.OlMo (see Sec. 0, the mass of the formed black 
hole should be approximately ^ 0.99Afo- Assuming that 
the area of the apparent horizon is equal to that of the 
event horizon, the non-dimensional spin parameter of the 
black hole is defined hy q = Jbh/-^^bh where Jbh and 
A/bh are the angular momentum and the ADM mass of 
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the black hole. Then the irreducible mass is defined in 
terms of q and Mbh by 



/l + fl - 0^)1/2 

A/to = y 2 ^^^^ 

Equation (EH) imphes that for Mirr/Mo = 0.9-0.97, 
g = 0.8-0.4. This indicates that the small error in the 
estimation for Mi„ leads to a large error in q estimated 
from the area of the apparent horizon. To reduce the 
error, the simulation should be performed with a better 
grid resolution. For this purpose, an adaptive mesh re- 
finement technique will be helpful |6l| , but the simulation 
with such technique is beyond the scope of this paper. 

The value of q is also approximately determined from 
the following manner. As shown in Sec. the angular 
momentum is dissipated by ^ 15% by gravitational radi- 
ation, while the ADM mass decreases by ~ 1% through- 
out the simulation. As listed in Table II, the initial value 
of g is « 0.9. Therefore, the value of q in the final stage 
should be « 0.75 for which Mi„ ~ 0.91 Mbh- Assuming 
that this is the correct value, the error in the irreducible 
mass of the black hole determined from the area of the 
apparent horizon is within 5%. 

Given the spin parameter q — 0.7-0.8 and the ADM 
mass A/bh ~ 2.8Mq, the frequency of the fundamental 
quasinormal mode is about 6.5-7(2.8M0/Mbh) kHz [6^. 
We will show the ring-down gravitational waves with this 
frequency in Sec. 



D. Implication to SGRBs 

As shown in sections IIVBI and IIV CI (also illustrated 
in previous papers pll l57|'). merger of binary neutron 
stars can produce a system composed of a Kerr black 
hole and a hot torus of mass ^ O.OIMq if certain condi- 
tion is satisfied. To summarize, there are following two 
scenarios to achieve the formation of a massive disk with 
Md > O.OIM0: (i) In the first scenario, a hypermassive 
neutron star is required to be formed first. Then, by some 
mechanisms, the angular momentum is transported from 
the inner region to the outer envelop which subsequently 
forms a hot disk (torus) , while the hypermassive neutron 
star eventually collapses to a black hole either by angular 
momentum dissipation due to gravitational radiation or 
by angular momentum transport probably due to mag- 
netic effects |23|. In this scenario, the estimated torus 
mass will be ^ O.OI-O.OSMq in the absence of strong 
magnetic effects, while in its presence, the mass can be 

0.05Mo |5a|57|. (u) In the second scenario, a black 
hole is promptly formed after the onset of merger in an 
unequal-mass binary of sufficiently small mass ratio Qm- 
Then, a torus is formed from the less massive neutron 
star which is tidally elongated at the merger and subse- 
quently constitutes a hot and geometrically thick torus. 
In this case, the torus mass depends strongly on the value 
of Qm and nuclear EOSs. 



In |57l | , a scenario through formation of a hypermassive 
neutron star with strong magnetic fields is described. In 
this case, a large accretion rate from the torus to the 
black hole with M > 5Mq/s is expected. In the rest of 
this section, we focus mainly on other scenarios with no 
strong magnetic fields; in particular, we focus on scenario 

(ii)- 

In the following, we assume that the disk mass is 
Md ~ O.OI-O.IM0 and mass accretion rate is Af ~ 0.1- 
1Mq/s. These are plausible values for a central engine 
of SGRBs since the life time of torus becomes tdur ~ 10- 
1,000 ms with such choice. The transport time of the 
angular momentum used here for M is assumed to be 
determined by an appropriate value of the viscosity or 
moderately strong magnetic fields as in the a-viscosity 
model (63.] . 

The gravitational binding energy W of disk with mass 
~ 0.01-0. lAf0 in the black hole spacetime is 

^ ^ AfBHAfd 
Td 

„ ^ .1 /IOA/bhW Md \ 
« 2 X 10 ergs , (32) 

where denotes an averaged disk radius. Due to vis- 
cous dissipation in the accretion disk, the kinetic energy 
is converted to the thermal energy. Equation l|32|l in- 
dicates that the total thermal energy produced during 
the accretion will be ~ 10^^-10^^ ergs, which is lO^-lO'^ 
times larger than the total radiated energy in SGRBs 
(for the value after correction of a beaming factor) 0]. 
Therefore, a relativistic fireball of sufficient energy may 
be formed if there is a mechanism by which a fraction 
of the thermal energy is converted to the fireball. In the 
past several years, scenarios f or p roducing such fireball 
have been indeed proposed |Tol lea IgI lei IgI IgtL IM l69l| . 
We describe a scenario following a semianalytic calcula- 
tion je^ • More detailed numerical simulations have been 
performed recently Ib^ . According to these numeri- 
cal works, the neutrino luminosity can be by a factor of 
^ 10 larger than that in 67] if the mass of disk is suffi- 
ciently large ^ O.O5Af0. Thus, the luminosity described 
below may be considered as a conservative value. 

Because of its high temperature (typically ~ 10^^ K) 
and density (- 10" g/cm^; cf. Figs □ |1 and Eld)- 
(f)), the torus radiates strongly in thermal neutrinos. 
The opacity inside the torus (considering only neutrino 
absorption and scattering interactions with nucleons) is 
K - 7 X 10~i^(r/10" K)2 cm2 g-i [el 113 where T 
is the temperature. To estimate the optical depth, we 
define the surface density of torus by 

^i^^y) ^ / pu^^f^dz, (33) 

where the integral is carried out along lines of a; = 
y ^constant. For Md ^ O.O1A^0 (e.g., for models 
APR1317, APR127175, APR1218, SLyl315, SLyl25155, 
SLyl216), S ^ 10^"^ g/cm^, so that the optical depth 
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of the neutrinos ~ kS is ^ 1 for r <, 20 km « 5Mbh- 
This optical depth is in approximate agreement with the 
model presented in (Note that in a high viscosity 

case, the significant increase of the thermal pressure by 
the viscous dissipation may enforce the torus to expand 
and to decrease the optical depth below unity 

In the optically thick case, approximate neutrino lu- 
minosity may be estimated in the diffusion limit 70] as 
Li, ~ "^r^F where F is the neutrino flux from the torus 
surface, approximately estimated by 



7N^ aT'^ 
3 kS 



(34) 



Here, a denotes the Stefan-Boltzmann constant and A^^ 
is the number of neutrino species, taken as 3. Then, the 
neutrino luminosity is expressed as 



~ 2 X 10^^ ergs/s 



10 km 



10" KJ VlO" g/cm^ 



, (35) 



which is only slightly smaller than the neutrino Edding- 
ton luminosity Hence, our numerical results suggest 
formation of a hot, hyperaccreting torus which is opti- 
cally thick to neutrinos. A model for the neutrino emis- 
sion in a similar flow environment with comparable Ly is 
provided in (67j as a neutrino-dominated accretion flow 
(NDAF). 

According to the model of 67], the luminosity due to 
neutrino-antineutrino {vv) annihilation is [6^ 



L^i, - 10^^-10^" ergs/s, 



(36) 



for M - O.l-lMg/s. (The work in j6J,|63 also indicates 
that for M < 0.1Mq/s, the luminosity steeply decreases 
far below Li^p = 10^^ ergs/s.) In our model, the black 
hole is rapidly rotating with q ~ 0.75, and hence, the 
luminosity may be enhanced by a factor of several due 
to a GR effect [g^]- Also, in a more careful estimation 
for the neutrino annihilation and for the neutrino opacity 
than those in 67], the value of could be increased by 
a factor of - 10 l6i(7l|. 

Because of the thick geometry of the torus, pair fD 
annihilation will be most efficient near the z-axis and just 
above the surface of the inner region of the torus |68L uM 
for which the density is much smaller than that of the 
torus. Consequently, the baryon- loading problem will 
be escaped, and hence, a strongly relativistic fireball is 
likely to be produced. 

Aloy et al. [t^ simulate the propagation of jets pow- 
ered by energy input along the rotation axis (as would be 
supplied by the annihilation). They find that if the 
half-opening angle of the energy injection region is mod- 
erately small (^ 45°) and the baryon density around the 
black hole is sufficiently low, jets with the Lorentz factors 
in the hundreds can be produced given an energy input 
i^p ~ 10'*^-10^° ergs/s lasting ~ 100 ms. In such case. 



the conversion rate from the pair annihilation luminosity 
to the jet energy is a few x 10% in their results. They also 
indicate that the duration of SGRBs may be ~ 10 times 
longer than the duration of the energy input because of 
the differing propagation speeds of the jet head and tail; 
namely, duration of the energy supply with ~ 10 ms may 
be sufficient for explaining a SGRB with duration ~ 100 
ms. 

Our numerical results, along with the accretion fiow 
(NDAF) [s^ and jet propagation models of 72], thus 
suggest that the system composed of a black hole and hot 
torus with Md ^ O.OIM© and accretion rate M ~ 0.1- 
IMq/s is a possible candidate for the central engine of 
SGRBs. Since the lifetime of the torus is ^ 10-100 ms 
in our model, the total energy of the annihilation is 
E^i, - 10"** ergs(Md/0.01MQ). Typical total energy of 
SGRBs are several xlO'** ergs 0, and hence, a system 
of the disk mass ^ several XO.OIM0 is likely to be ap- 
propriate for powering SGRBs as long as the emission is 
beamed. Probably, beaming is encouraged by the thick 
geometrical structure of the torus [t^ • 

The above discussion suggests that if the mass ratio is 
sufficiently small, prompt formation of a black hole may 
be a possible scenario for producing the central engine 
of SGRBs. In particular, for a small value of Qi\i <, 0.8 
with the SLy EOS, a large SGRB energy is expected. 
Unfortunately, binary neutron stars of small mass ratio 
with Qm < 0.9 have not been found so far Hence, 
it is not clear at present whether the formation rate of 
such unequal-mass binary is large enough for explaining 
the event rate of SGRBs. This is a weak point in this 
scenario. 

Through formation of a hypermassive neutron star 
with moderately strong magnetic fields, a system com- 
posed of a black hole and a disk of mass ~ O.O1-O.O3M0 
may be formed as discussed in Sec. IIV Bl The estimation 
described in this section is also applied for such system, 
and indicates that a central engine of SGRBs of rela- 
tively low total energy ^ several x 10''^-10'*^ ergs could 
be formed. If the time scale of the angular momentum 
transport by the magnetic effects (or viscosity) is longer 
than emission time scale of gravitational waves, a black 
hole is formed (i.e., a SGRB should be generated) after 
a longterm emission of quasiperiodic gravitational waves 
from the hypermassive neutron star. As we will discuss 
in Sec. such gravitational waves may be detected by 
advanced laser-interferometric detectors. If quasiperiodic 
gravitational waves and subsequent SGRB are detected 
coincidently in the same direction with a small time lag, 
thus, this scenario may be confirmed. 

The results for the disk mass in this paper suggest 
that a merger between a black hole and a neutron star 
could form a system composing a black hole and sur- 
rounding massive disk by tidal disruption 73]. In- 
deed, if the mass of the black hole is small enough 
[Mbh ^ 4M0(i?o/lO km)3/2 for g = and Mbh ^ 
20Mo(i?o/10 km)3/2 for q = 0.9 where q and Rq are 
a spin parameter of the black hole and a neutron star 
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radius ^7^], the neutron star will be tidally disrupted 
before reaching the ISCO. If the event rate (though it is 
not clear at present because black hole-neutron star bi- 
nary has not been observed so far) is large enough, such a 
merger could be a promising source for producing a cen- 
tral engine of SGRBs. However, a detailed simulation in 
full general relativity is necessary to confirm this scenario 
since the disk mass depends crucially on the location of 
the ISCO. 



V. GRAVITATIONAL WAVEFORMS 

A. Gravitational waves from hypermassive neutron 

stars 

1. Waveform and luminosity 

In Fig. 1141 we present gravitational waveforms in 
the formation of hypermassive neutron stars (models 
APR1313 and APR1414) as a function of a retarded time. 
Throughout this paper, the retarded time tret is defined 
by t — 7'obs where robs is the coordinate radius of the 
wave extraction. In the early phase (tj-ct ^ 2 ms), grav- 
itational waves associated with the inspiral motion are 
emitted, while for ij-ct 2 ms, those by the rotating and 
oscillating hypermassive neutron star are emitted. In the 
following, we focus only on the waveforms for t^ct <; 2 ms. 

For model APR1313, a hypermassive neutron star of 
ellipsoidal shape is formed after the merger sets in. As a 
result, quasiperiodic gravitational waves with an approx- 
imately constant frequency « 3.2 kHz are emitted. Also, 
the amplitude remains approximately constant for t ^ 5 
ms. These properties are essentially the same as those 
found in 21] for models SLyl313 and SLyl25135. 

For model APR1414, on the other hand, the ampli- 
tude decreases with time in particular for < > 6 ms. 
This reflects the fact that the ellipticity of the formed 
hypermassive neutron star steeply decreases for i ^ 6 ms. 
This indicates that only less massive binaries can produce 
a longterm emitter of quasiperiodic gravitational waves. 
For model APR1414, the frequency is not constant ei- 
ther, but modulates with time as in model SLyl35135 
[2l|. The reason is that the formed hypermassive neu- 
tron star quasiradially oscillates with a large amplitude 
and the characteristic frequency varies with the change 
of the characteristic radius. Due to this, the shape of 
the Fourier spectra for models APR1313 and APR1414 
are different even qualitatively (cf. Fig. I17|l . Namely, 
for M <i Mthr, a slight difference in the total ADM mass 
results in a significant difference in gravitational wave- 
forms. 

In Fig. El the emission rates of the energy and the an- 
gular momentum by gravitational radiation are shown for 
models APR1313 (solid curves) and APR1414 (dashed 
curves). In the inspiral phase for t^ct ^ 2 ms, the emis- 
sion rates increase with time (besides initial unphysical 
bump associated with the conformal flat initial condition 



in which gravitational waves are neglected) , since the am- 
plitude and the frequency of the chirp signal increase. 
Then the peak is reached at t^ct ~ 2 ms. The height 
of the peak is larger for model APR1414 since the com- 
pactness of each star is larger. After the peak is reached, 
the emission rates once quickly decrease since the merged 
object becomes a spheroidal transient object. However, 
because of its large angular momentum and stiff EOS, 
the formed hypermassive neutron star soon changes to a 
highly ellipsoidal object which emits gravitational waves 
of a large amplitude. The luminosity from the ellipsoidal 
neutron star is as high as the first peak at tj-d ~ 2.5 
ms for model APR1414 and at t^ct ~ 3.5 ms for model 
APR1313. 

After the second peak, the emission rates of the en- 
ergy and the angular momentum via gravitational waves 
gradually decrease with time, since the degree of the non- 
axial symmetry decreases. However, for model APR1313, 
the decrease rates are small and the emission rates at 
^rot ~ 10 ms remain to be as high as that in the late 
inspiral phase as dE/dt ~ 6 x 10'^'* erg/s and dJ/dt ~ 
6 X 10^" g cm^/s^ reflecting the high ellipticity of the 
hypermassive neutron star. The angular momentum at 
i ~ 10 ms is J ~ 0.7 Jo ~ 4 x 10^^ g cm^/s. Assuming 
that the emission rate of the angular momentum remains 
~ 5 X 10^° g cm^/s and that emission stops when about 
half of J is dissipated, the emission time scale is approx- 
imately evaluated as J/ {2dJ/dt) ^ 40 ms. If the nonax- 
isymmetric structure is maintained and dJ/dt does not 
vary much, thus, the hypermassive neutron star will col- 
lapse to a black hole within ~ 50 ms even in the absence 
of other dissipation processes (cf. the discussion in the 
last paragraph of IIV B l")) . 

For model APR1414, on the other hand, the emission 
rates decrease quickly since the hypermassive neutron 
star relaxes to an approximately axisymmetric spheroid 
for t k;, 6 ms. For this model, the dissipation time scale 
of the angular momentum is much longer than 50 ms at 
t = 10 ms. Therefore, other dissipation processes such as 
magnetically induced angular momentum transport will 
trigger the collapse to a black hole. 

By the time integral of dE/dt and dJ/dt, the total 
energy and angular momentum radiated are computed 
and found to be about 0.03A/o and 0.30Jo for model 
APR1313 and 0.03Mo and 0.26Jo for model APR1414, 
respectively. This indicates that the angular momentum 
is significantly dissipated, illustrating that the angular 
momentum dissipation plays an important role in the 
evolution of the system. To confirm that the radiation 
reaction is followed in the simulation, we display J{t) and 
Jo — A(i) as a function of time for models APR1313 and 
APR1414 in Fig. QSl This shows that the angular mo- 
mentum computed from Eq. (|20|l agrees approximately 
with Jq — AJ (within ~ 2% error), proving that radiation 
reaction is computed with a good accuracy. 
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FIG. 14: Gravitational waveforms, R+ and Rx , (a) for model APR1313 at robs = 36Afo and (b) for model APR1414 at 
robs = 32Mo. 




FIG. 15: Energy and angular momentum emission rates, 
dE/dt and dj/dt, of gravitational waves for models APR1313 
(solid curves) and APR1414 (dashed curves). The units are 
erg/s and g cm^/s^, respectively. 
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FIG. 16: The solid curve denote evolution of angular momen- 
tum for models (a) APR1313 and (b) APR1414. The dashed 
curves denote Jo — A J where Jo denotes the angular momen- 
tum at t = 0. 



2. Fourier spectrum 

In the real data analysis of gravitational waves, a 
matched filtering technique is employed. In this 
method, the signal of the identical frequency can be ac- 
cumulated using appropriate templates. As a result, the 
effective amplitude increases by a factor of A^^/^ where 
N denotes an approximate number of the cycle of gravi- 
tational waves for a given frequency. 

To determine the characteristic frequency of gravita- 
tional waves, we carry out a Fourier analysis. In Fig. 
1171 the power spectrum dE/df is presented for mod- 
els APR1313 and APR1414. Since the simulations were 



started with the initial condition of the orbital period 
~ 2 ms (i.e., frequency of gravitational waves is ~ 1 
kHz), the spectrum of inspiraling binary neutron stars 
for / < 1 kHz cannot be correctly computed. Thus, only 
the spectrum for / 1 kHz should be paid attention. As 
a plausible spectrum for / ^ 1 kHz, we plot the Fourier 
power spectrum of two point particles in circular orbits 
in the second post Newtonian approximation (the dotted 
curve) (the third post Newtonian terms does not sig- 
nificantly modify the spectrum since their magnitude is 
~ 0.01 of the leading-order term). 

Figure El shows that a sharp characteristic peak is 
present at / « 3.2 and 3.8 kHz for models APR1313 and 
APR1414, respectively. This is associated with quasiperi- 
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FIG. 17: Fourier power spectrum of gravitational waves 
dE/df for models APR1313 (solid curve) and APR1414 
(dashed curve). Since the simulations are started when the 
characteristic frequency of gravitational waves is ~ 1 kHz, 
the spectrum for / < 1 kHz cannot be presented. The dotted 
curve in the panel denotes the analytical result of dE/df in 
the second post Newtonian and point-particle approximation 
by which the real spectrum for f ^ 1 kHz is approximated. 

odic gravitational waves emitted by the formed hyper- 
massive neutron stars. Two side-band peaks are present 
at / « 3.2 and 4.7 kHz for model APR1414. Thus, 
the spectral shape is qualitatively different from that for 
model APR1313. The reason is that the amplitude of the 
quasiradial oscillation of the hypermassive neutron star 
is outstanding and the characteristic radius varies for a 
wide range for model APR1414, inducing the modulation 
of the wave frequency. 

An effective amplitude of gravitational waves observed 
from the most optimistic direction (which is parallel to 
the axis of the angular momentum) is proportional to 
yj dE/df in the manner 

/leff = Y^|i?+|2 + |i?x|2/ 

= i.8xio-f^;^V'7H^V) 

VlO^i erg/Hz / V ^ / 

where r denotes the distance from the source, and i?+,x 
are the Fourier spectrum of R+^y.- In Fig. we show 
/loff as a function of / for a hypothetical distance of 50 
Mpc. This shows that the effective amplitude of the peak 
is ~ 3 times larger than that at ^ 1.3-1.5 kHz which 
corresponds to the frequency of the last inspiral motion. 

For model APR1313, furthermore, the amplitude of 
the peak in reality should be larger than that presented 
here, since we stopped simulations at t ^ 10 ms to save 
the computational time, and hence, the integration time 
(^ 10 ms) is much shorter than the realistic value. Ex- 
trapolating the decrease rate of the angular momentum, 
the hypermassive neutron star will dissipate sufhcient an- 
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FIG. 18: Non-dimensional effective amplitude of gravitational 
waves from hypermassive neutron stars for models APR1313 
(solid curve) and APR1414 (dashed curve). The assumed 
distance is 50 Mpc. The dotted line denotes the planned 
noise level of the advanced LIGO. 



gular momentum by gravitational radiation until a black 
hole or a spheroidal star is formed. As indicated in Sec. 
IV A II the duration of the angular momentum dissipation 
would be 50 ms. Thus, we may expect that the emis- 
sion will continue for such time scale and the effective 
amplitude of the peak would be in reality amplified by a 
factor of ~ 5^/^ « 2 to be ~ 6 x 10"^^ at a distance of 
50 Mpc. Although the sensitivity of laser-interferometric 
gravitational wave detectors for / > 1 kHz is limited by 
the shot noise of the laser, this value is by a factor of 
3-4 larger than the planned noise level of the advanced 
LIGO « 10-2i-5(//l kHz)3/2 13 . It wiU be interesting 
to search for such quasiperiodic signal of high frequency 
if the chirp signal of gravitational waves from inspiral- 
ing binary neutron stars of distance r ^ 50 Mpc are 
detected in the near future. As discussed in [23|, the de- 
tection of such quasiperiodic signal leads to constraining 
nuclear EOSs for neutron star matter. Furthermore, as 
mentioned in Sec. IIVDI the detection may also lead to 
confirming a scenario for producing the central engine of 
SGRBs. 

On the other hand, for model APR1414, the quasiperi- 
odic gravitational waves damp at t ~ 6 ms, and hence, 
the expected amplitude is hcs ^ 2 x 10~^^ at a distance 
of 50 Mpc with f ^ 4: kHz. If the required signal to noise 
ratio of gravitational waves for confirming the detection 
is '--^ 5 in the advanced LIGO, such signal will be detected 
only for an event within the distance ~ 10 Mpc. Since 
the predicted event rate of the merger for such distance is 
<i 0.01 per years frequent detection is not expected. 
However, the fortunate detection may lead to confirming 
a scenario for producing the central engine of SGRBs. 



24 




-0.5 



solid curve: r = 27M 
dotted curve: fitting 
J I I I I I I 



2.2 



2.4 



2.6 



t-r 



obs 



ms, 



X 




-0.5 



2.8 



solid curve: r = S9M 
dotted curve: fitting 

J I I I I I \ I 1 



3.2 



t-r 



obs 



ms, 



3.4 



FIG. 19: The same as Fig. [H but for models APR135165 (left) at robs = 27Afo and model SLy 125155 at robs = 29Mo (right). 
The lower panel is the enlargement of the upper panel for ring-down gravitational waves associated with a quasinormal mode 
of the formed black hole. The formation time of the apparent horizon tAH is 2.38 ms for model APR135165 and 3.04 ms for 
model SLyl25155, and hence, the quasinormal mode is excited for t — robs ^ tAH. The thick dotted (blue) curves denote the 
fitting formulae described by Eq. 13811 . 



B. Gravitational waves in the black hole formation 

In the prompt formation of a black hole, ring-down 
gravitational waveforms associated with a quasinormal 
mode is emitted. Figure E| shows i?+ for model 
APR135165 and Rx for model SLyl25155 as a function of 
retarded time. In the early time (tret ^ 2 ms for model 
APR135165 and U-^t 2.3 ms for model SLyl25155), 
gravitational waveforms are determined by the inspiral- 
ing motion. After this, a merger waveform, which is emit- 
ted due to hydrodynamic interaction between two stars, 
is seen for a short time scale ^ 0.5 ms. The merger wave- 
form is emitted for a longer time for model SLyl25155 
simply because the time from the onset of the merger to 
formation of a black hole is longer. Finally, a black hole 
is formed and the waveforms are determined by the fun- 
damental quasinormal mode. The formation time of the 
apparent horizon tAH is 2.38 ms for model APR135165 
and 3.04 ms for model SLyl25155. Assuming that the 
formation time of the event horizon is approximately the 
same, ring-down waveforms associated with the quasi- 
normal mode should be induced only for tret ^ ^ah- In 
our numerical results, this condition indeed holds, and 
therefore, we have confirmed that the quasinormal mode 
is extracted [76j . 

For models APR1416, APR135165, APR1317, and 
SLyl25155, the initial values of the non-dimensional an- 
gular momentum parameter qq is « 0.9. The energy and 
the angular momentum are dissipated by gravitational 
waves by ~ O.OIA/q and O.ISJq until formation of the 
black holes (cf. Fig. I20|l . Taking into account the fact 
that a tiny fraction of mass and angular momentum is 



distributed to the surrounding disk, the spin parameter 
of the formed black hole should be ~ 0.75. The dotted 
curve in Fig. 1191 denotes a model of ring-down gravita- 
tional waveforms 

^g-(t-tAH)/t. cos(2^/q„^t + (5), (38) 

where A denotes the maximum amplitude and S a phase 
constant, /qnm and td are the frequency and damping 
rate of ring-down gravitational waves oil = m = 2 mode, 
which are approximately given by [g^ ] 

/qnm - 10.8 (^^^ [1 - 0.63(1 - qf-^] kHz, (39) 



^/qnm 

For M = 2.9Af0 and q = 0.75, /qnm = 6.5 kHz and td = 
0.183 ms which are used for the fitting to the waveform 
of model APR135165. For model SLyl25155, we plot 
the curve with M — 2.7Mq and q = 0.7, and hence, 
/qnm — 6.7 kHz and td = 0.163 ms. This illustrates that 
with M = 2.7-2.9Af0, the frequency will be in a small 
range between 6.5 and 7 kHz. 

In Fig. ^1 the fitting formulae are plotted together. 
It shows that the computed waveforms agree approxi- 
mately with the fitting curves. This implies that gravita- 
tional waves numerically extracted carry the information 
of the formed black hole correctly. The maximum value 
of R+ X in the ring-down phase is ~ 0.15 km for model 
APR135165 and - 0.25 km for model SLyl25155, im- 
plying that the maximum amplitude is ~ 1-2 x 10~^^ at 
a distance of 50 Mpc [cf. Eq. <(T7|l ]. These values de- 
pend weakly on the mass ratio. The amplitude in model 
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SLyl25155 is slightly larger than that in the APR EOS. 
This is probably because the accretion rate of the matter 
which can excite the quasinormal mode is larger in the 
SLy EOS. 

Since the frequency is too high and the amplitude is 
too small, it will be difficult to detect ring-down gravita- 
tional waves even by advanced laser-interferometric grav- 
itational wave detector. Only in the case that a merger 
happens in our local group, it may be detected. However, 
the expected event rate is less than 0.001 /yrs 12]. 
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FIG. 20: Evolution of the total radiated energy /S.E and an- 
gular momentum A J for models APR1416 (solid curves), 
APR135165 (dashed curves), and APR1317 (long dashed 
curves) 

Figure 1^ shows the evolution of radiated energy and 
angular momentum by gravitational waves for models 
APR1416, APR135165, and APR1317. R is found that 
~ 1% of the total mass energy and ~ 15% of the angular 
momentum are dissipated by gravitational waves. These 
gravitational waves are mostly emitted in the last one 
inspiral orbits and the contribution from the ring-down 
phase is negligible. (For model SLyl25155, the results 
are quantitatively the same as those of APR models.) 
The significance of the dissipation of the angular mo- 
mentum in the last orbit should be emphasized since it 
affects the mass of disks surrounding the formed black 
hole. Another important property found from Fig. 1201 
is that for smaller mass ratios, the radiated energy and 
angular momentum are smaller. This results from the 
fact that with the smaller mass ratios, the tidal elonga- 
tion sets in at a larger orbital separation, and hence, the 
inspiral waveform shuts off earlier 0, |^ . 

VI. SUMMARY 

As an extension of our previous work plj . we per- 
formed fully general relativistic simulations for the 
merger of binary neutron stars adopting stiff nuclear 



EOSs. We focus particularly on the formation of a black 
hole and surrounding disk in this paper. The following is 
the summary of the new results described in this paper. 
The summary about the outcome after the merger is also 
described in Fig. [2] 

1. Threshold mass for black hole formation 

In the APR EOS which is mainly adopted in this pa- 
per, the threshold mass for prompt formation of a black 
hole is Mthr ~ 2.8-2.9Mo which is about 30-35 % larger 
than the maximum mass for spherical neutron stars with 
identical EOS, Mgph- Collecting the results obtained in 
|2ll |. we conclude that the following relation holds for 
stiff nuclear EOSs such as EPS, SLy, and APR EOSs: 
Afthr = 1.3-1.35Msph. 

2. Formation of hypermassive neutron star (M < ^fthr) 

If M < Mthr, a hypermassive neutron star is formed 
as an outcome of the merger. If the mass is not close to 
Afthri the outcome has an ellipsoidal shape, and hence, 
is a strong emitter of quasiperiodic gravitational waves. 
Because of angular momentum dissipation by gravita- 
tional waves, the ellipsoidal hypermassive neutron star 
may collapse to a black hole within ^ 50 ms, as found in 
[2l|, even in the absence of any other transport mecha- 
nism. On the other hand, if the ellipticity reduces to zero 
within ~ 50 ms, the hypermassive neutron star will settle 
to a spheroid in a stationary state. Then, the collapse to 
a black hole will be triggered by other mechanisms such 
as the magnetic braking and the MRI [ssl . Is^ Is^ . 

Since a torque works from the ellipsoidal hypermassive 
neutron star, matter in the outer envelop receives angular 
momentum from the central part. This process helps 
producing disk, which may be an accretion disk around 
a black hole eventually formed. If gravitational radiation 
triggers the collapse, the expected disk mass is ^ 0.01- 
O.O3Af0 for merger of equal-mass as well as unequal-mass 
binaries for Qm ^ 0.85 with the black hole mass ~ 2.6- 
2.8M0 and spin q < 0.7. Such system may be a central 
engine of SGRBs of relatively small burst energy |^ . 

Alternatively, collapse to a black hole may be triggered 
by other mechanisms such as the magnetic braking and 
the MRI ^55, 56, 57] . In the collapse by these magnetic 
effects, the mass of a torus surrounding the black hole will 
be larger than O.O5M0 because an efficient angular mo- 
mentum transport works, as illustrated in ^56, 57]. Fur- 
thermore, the MRI generates shocks in the torus, heating 
up the material to ~ 10^^ K js^]. Therefore, this sce- 
nario is more favorable for producing a central engine of 
SGRBs. 

If the time scale of the angular momentum transport in 
the hypermassive neutron stars is longer than ^ 50 ms, 
a black hole or a spheroidal hypermassive neutron star is 
formed after a longterm emission of quasiperiodic gravi- 
tational waves. Such gravitational waves may be detected 
by advance laser-interferometric detectors (see summary 
6) . If a black hole is subsequently formed, it may also pro- 
duce a SGRB. Thus, if quasiperiodic gravitational waves 
and subsequent SGRB are detected coincidently in the 
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FIG. 21: Summary about the outcome after the merger. "HMNS", "GW", "B-field", and "J-transport" are hypermassive 
neutron star, gravitational wave emission, magnetic field, and angular momentum transport, respectively, "small disk" , "massive 
disk", and "heavy disk" imply that the disk mass is Md < O.OIM©, O.OIM© ^ Ma <, O.OSMq, and Md ^ O.OSM©, respectively 
The solid arrow denotes the route found in this paper and the dashed arrow is the possible route based on the results of pTII 
and the speculation. There are three possible fates of the ellipsoidal hypermassive neutron star depending on the EOS and the 
time scale of dissipation and angular momentum transport processes. We note that a spheroid is formed probably only for the 
APR EOS. See the summary 1-4 for details. 



same direction with small time lag, the merger scenario 
through formation of a hypermassive neutron star for the 
central engine of SGRBs may be confirmed. (However, 
for detecting gravitational waves, the event should be 
within ^ 50 Mpc; see summary 6). 

In the presence of strong magnetic fields with J> 10^^ 
G, a black hole may be formed before a longterm emission 
of quasiperiodic gravitational waves from the hypermas- 
sive neutron star 56, 57]. In this case, a smaller value of 
/icff for the quasiperiodic gravitational waves is expected, 
and hence, a smaller-distance event is required for confir- 
mation of this scenario by gravitational wave detection. 

3. Formation of hypermassive neutron star (M ^ Mthr) 

For M only slightly smaller than Mthr, the hypermas- 
sive neutron star relaxes to a spheroidal state in ~ 10 ms 
since its small angular momentum and compact structure 
lead to a too small ratio of the rotational kinetic energy 
to the gravitational binding energy to form an ellipsoid 
|58l | . Here, the small angular momentum is a result from 
the effects that the angular momentum is significantly 
dissipated by gravitational waves soon after the onset of 
the merger and that the angular momentum transport 
works efficiently during an early quasiradial oscillation 
with a large amplitude in the hypermassive neutron star. 
Since the emissivity of gravitational waves by the hy- 
permassive neutron stars of small ellipticity is not high, 
it may survive for more than 50 ms in contrast to the 
smaller-mass case. Collapse to a black hole will be trig- 
gered by other mechanism such as the magnetic braking 
and the MRI 55, 56, 5'^. After the collapse, disk of mass 
^ 0.05A/o could be formed H^H^, and hence, such sys- 
tem is also a candidate for the central engine of SGRBs. 



In this scenario, the amplitude of quasiperiodic gravita- 
tional waves is not as high as that for the smaller-mass 
case, and hence, frequent detection of this signal is not 
expected (see summary 6). 

The formation of a spheroid is found only for the APR 
EOS (in IP in which the SLy and FPS EOSs are used, 
we have not found the formation of spheroids). The likely 
reason in this difference is that the APR EOS is stiffer 
than others. Namely, the hypermassive neutron stars can 
be more compact with this EOS [cf. Fig. Olb)] escaping 
the collapse to a black hole. In such compact state, the 
gravitational binding energy could be large enough to 
reduce the ratio of the rotational kinetic energy to the 
binding energy below the threshold value of the formation 
of an ellipsoid. 

4. Prompt formation of black hole {M > A/thr) 

In the black hole formation from the nearly equal-mass 
merger, most of mass elements are swallowed into the 
horizon, and hence, the disk mass around the black hole 
is much smaller than O.OIMq. However, the disk mass 
steeply increases with decreasing the value of Qm ; we find 
an empirical relation in which the disk mass is written 
approximately by Mdo + Mdi{l — QmY with p = 3-4 for 
given EOS and mass for a range 0.7 <^ Qm < 1- 

For Qm ^ 1, the disk not only has small mass but 
also is geometrically thin. Therefore, the merger will not 
produce a central engine of SGRBs for nearly equal-mass 
case. On the other hand, for a sufficiently small value 
of Qm {Qm ^ 0.75 for the APR EOS and Qm ;^ 0.85 
for the SLy EOS), the disk mass is larger than O.OIM©. 
Furthermore, the disk is geometrically thick and the ther- 
mal energy is large enough (typically ~ 1-2 x 10^^ K) 



27 



for producing the large amount of thermal neutrinos. In 
addition, the black hole spin is large q ~ 0.75. There- 
fore, prompt formation of a black hole in the merger of 
unequal-mass binary neutron stars with a small value of 
Qm is a scenario for formation of a central engine of 
SGRBs. Unfortunately, binary neutron stars of small 
mass ratio with Qm < 0.9 have not been found so far 
|49|. Hence, it is not clear at present whether the merger 
rate of such unequal-mass binary is large enough for ex- 
plaining the event rate of SGRBs. The present results 
also indicate that the merger between a neutron star and 
a black hole with small mass will be a possible candidate 
for producing central engine of SGRBs. 

5. Gravitational waves from black hole 

The non-dimensional angular momentum parameter 
{J /A'P) of the formed Kerr black hole is ^ 0.75 in the 
prompt formation case. Then, for the system of mass 2.7- 
2.9M0, the frequency of ring-down gravitational waves 
associated with the fundamental quasinormal mode of 
Z = TO = 2 is ~ 6.5-7 kHz. We extract the ring-down 
gravitational waveform and confirm this frequency. The 
amplitude of gravitational waves is ^ 1-2 x 10~^^ at a 
distance of 50 Mpc which is too small to be detected 
even by advanced laser-interferometric detectors unless 
the merger event happens within the local group of galax- 
ies. 

6. Gravitational waves from hypermassive neutron star 
The effective amplitude of quasiperiodic gravitational 



waves from hypermassive neutron stars of ellipsoidal 
shape can be ^ 5 x 10~^^ at a distance of 50 Mpc 
with the frequency 3-3.5 kHz. This property agrees 
with that found in [2lll2^ . and hence, we conclude that 
this fact holds irrespective of stiff nuclear EOSs with 
Ms-ph ^ 2 A/0. Quasiperiodic gravitational waves may 
be detected by advanced laser-interferometric detectors. 
Detection will lead to constraining the nuclear EOS as 
discussed in |2^ . On the other hand, for M only slightly 
smaller than Mthr, the expected effective amplitude is 
/loff ~ 2 X 10"^"'^ at a distance of 50 Mpc with frequency 
/ ~ 4 kHz; only for an event with a distance within ^ 10 
Mpc, detection of such signal will be possible by the ad- 
vanced laser-interferometric detectors. 
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